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ABSTRACT 


Over  the  range  of  wavelengths  from  one  to  six  microns 
where  oxide  ceramics  are  partially  transparent,  radiant  energy 
emission  must  be  treated  as  a  volume  rather  than  a  surface 
phenomena.  As  a  result  the  absorption  coefficient  (mainly 
determined  by  composition)  and  the  scattering  coefficient 
(mainly  determined  by  micros true ture)  are  important  variables. 
In  addition,  temperature  gradients  affect  radiant  energy 
emission  so  that  the  effective  " emissivity"  under  non-iso- 
thermal  conditions  is  neither  a  material  constant  nor  a  sample 
constant. 

Analytical  relationships  have  been  derived  for  spectral 
radiant  energy  emission  as  a  function  of  absorption  coefficient, 
scattering  coefficient,  surface  temperature  and  temperature 
gradient.  Experimental  measurements  of  non- iso thermal  radiant 
energy  emission  of  AI2O3  and  SrTi(>3  samples  are  in  agreement 
with  these  calculations.  Analytical  relations  and  experimental 
measurements  for  diffuse  and  parallel  light  transmission  are 
reported.  Comparison  of  the  Hamaker  approximation  with  an 
exact  solution  show  that  the  approximation  is  suitable  for 
these  materials  and  conditions. 
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1.0  INTRODUCTION 


In  Part  I  of  the  present  study, ^ general  equations 
describing  emission,  transmission  and  reflectance  of  scattering 
and  absorbing  systems  were  derived  in  a  form  that  allows  deter¬ 
mination  of  the  necessary  material  constants  and  calculation 
of  emissivity  under  isothermal  and  non-isothermal  conditions. 
Methods  of  evaluating  conditions  under  which  these  equations 
are  applicable  were  derived,  and  a  simple  expression  for 
radiant  energy  transfer  when  conduction  in  the  solid  is  negli¬ 
gible  relative  to  the  radiant  energy  transfer  in  the  solid 
was  developed. 

Experimental  apparatus  was  developed  and  measurements 
made  of  isothermal  emissivity  of  several  refractory  oxides  and 
single  crystal  transmission  over  the  wavelengtn  range  1-15 
microns  at  temperatures  up  to  1350°C-  Emissivity  was  found  to 
depend  strongly  on  composition,  microstructure  and  wavelength. 
Different  grades  of  '•alumina*’  ceramics  range  from  about  0-2  to 
above  0-9  in  the  near  infra-red. 

Diffuse  transmission  measurements  permitted  calculations 
of  the  pertinent  material  constants  as  well  as  the  separate 
scattering  and  absorption  coefficients.  For  these  systems 
scattering  was  large  for  diffuse  radiation,  in  the  range  1000- 
5000  cm-1;  as  a  result  of  these  large  values  calculations  of 
emissivity  were  of  low  precision,  but  agreed  with  experimental 
measurements.  Results  showed  that  the  normal  emissivity  in 
these  systems  was  substantially  larger  than  the  hemispherical 
emissivity,  in  contrast  to  opaque  solids  where  they  are  nearly 
the  same . 

(2) 

In  Part  II  of  the  study  ,  emissivity  values  were  calcu¬ 
lated  for  specially  prepared  AI2O3,  MgO,  SiC>2  and  SrTiC>3  samples 
from  appropriate  reflectivities,  absorption  coefficients  and 
scattering  coefficients  derived  from  material  characteristics 
(index  of  refraction,  single  crystal  transmissivity,  pore  size 
and  concentration) . 

Experimental  measurements  of  emissivity  for  these  samples 
were  made  at  temperatures  up  to  120C°C  and  at  wavelengths  from 
one  to  fifteen  microns.  Calculated  and  measured  values  were 
in  good  agreement  for  AI2O3  and  SiC>2-  Absorption  coefficients 
of  the  MgO  and  SrTi03  samples  were  greater  than  the  single 
crystal  values,  leading  to  higher  measured  (than  calculated) 
emissivities .  Samples  of  commercial  alumina  had  higher 
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emissivities  than  the  specially  prepared  high  purity  samples. 

The  Baaaker  equations  for  non- iso thermal  emissivity,  which 
linearize  the  temperature  gradient,  were  compared  with  the 
exact  solution  by  numerical  integration.  Results  were  in  good 
agreement  for  all  cases  of  practical  interest  for  application 
to  ceramic  oxides.  Solutions  to  the  Hamaker  equations  were 
derived  for  new  boundary  conditions. 

Evaluation  of  the  rotating  blackbody  slot  system  used  for 
isothermal  emissivity  measurements  indicated  that  it  gave 
results  in  good  agreement  with  an  external  blackbody.  An 
apparatus  for  measuring  the  effective  emissivity  of  a  non- 
iso thermal  sample  was  constructed  and  tested. 

In  the  present  final  phase  of  the  study  our  major  effort 
has  been  devoted  to  further  evaluation  of  the  factors  affecting 
both  total  and  spectral  emissivity  from  ceramic  oxides  and 
carrying  out  experimental  measurements  of  non-isothermal 
emissivity  on  well  characterized  samples.  We  have  also 
extended  our  comparison  of  the  Ba maker  approximations  and  an 
exact  solution  to  actual  values  of  material  constants. 
Anticipating  our  results  we  can  say  that  we  have  experimentally 
observed  the  influence  of  non-isothermal  conditions  on  radiant 
energy  emission  from  ceramic  oxides  and  found  our  results  to 
be  in  reasonable  agreement  with  calculated  values. 

At  the  same  time  we  have  evaluated  analytically  the  trans¬ 
mission  of  radiant  energy  by  ceramic  oxides  containing 
scattering  centers,  and  made  experimental  measurements.  The 
low  radiant  energy  transmission  has  seriously  limited  our 
experimental  precision,  but  results  are  in  general  accord  vdth 
calculations . 

Overall,  we  now  have  a  sound  basis  for  estimating  the 
thermal  emissivity,  both  isothermal  and  non-isothermal,  from 
proper  material  characterization  and  for  evaluating  the  effect 
of  changing  material  characteristics  on  thermal  emissivity 
behavior . 
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2.0  NON- ISOTHERMAL  EMISSION  OF  RADIANT  ENERGY 


2.1  Total  Emission  of  Radiant  Energy 


We  shall  assume  that  a  sample  of  radiating  material 
has  a  specified  temperature  distribution  given  by 


T(x)  =  T  + 


.dT. 
(— )  x 
dx  o 


where  x  measures  the  distance  into  the  sample  from  the  surface. 
The  sample  has  an  absorption  coefficient  a,  scattering 
coefficient  S,  and  index  of  refraction  n.  Let  I*  (x,  6)  be 
the  specific  intensity  of  radiation  in  the  material,  where  9 
is  the  polar  angle  with  respect  to  the  +  x  axis  (Fig.  2.1) . 
Asstime  axial  symmetry. 


Vacuum  n  =  1 


Material  n  =  n 


I  (x>  0) 


+  x 


Figure  2.1  Diagram  of  Geometry  of  Radiation 

in  Absorbing  and  Scattering  Material 

The  equation  of  transfer  for  the  radiation  in  the  material 
is 

cos  9  =  -  (a  +  s)  I'(x,  9)  +  j  (x)  .  (2.1) 

j  (x)  is  the  emission  coefficient  for  radiation  in  the  material. 
It  has  two  terms,  one  due  to  thermal  emission  and  one  due  to 
scattering,  j  (x)  is  given  by 

j (x)  =  a  n2  o'T(x)4  +  S  I'  (x) ,  (2.2) 

where  a'  is  the  Stefan-Boltzmann  constant,  and 

TT 

I'(x)  =  2  tt  J  I'(x,  9)  sin  8  d0  (2.3) 

9=0 
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Now  assume  the  intensity  is  isotropic  over  each  hemis¬ 
phere  .  Take 


and 


I’  (x,  0)  *  I'  (x)  0  <  0  <  tt/2, 
I'  (x,  0)  *  I*  (x)  tt/2  <  9  <  tt. 


(2.4) 

(2.5) 


Then  we  easily  find  I'(x)  =  1/2  (I*  (x)  +1*  (x) ) .  The  flux 
is  given  by 

TT  _ 

ttF(x)  *  2ttJ  I  (x ,  0 )  cos  0  sin  0d  0  =  tt£i*  (x)  -  I'  (x)J 


0=  0 


In  fact,  in  terms  of  the  fluxes  I  and  J  we  have 


(2.6) 


I(x)  =  ttI*  (x)  and  J(x)  =  ttI*  (x)  . 

Now  multiply  the  equation  of  transfer  by  sin  0d6  and  in¬ 
tegrate  with  respect  to  0  from  0  to  tt-  Using  the  expressions 
(2)  and  (3)  for  I'(x,  0)  we  obtain,  after  some  simplification. 


=  -  4or  I '  (x)  +  4<*n2<j '  T  (x) 4 
dx 


(2.7) 


Similarly,  multiplying  the  equation  of  transfer  by 
cos  0  sin  8d0  and  integrating  as  before,  we  obtain  the  result 


4  di'  (x)  ,  .  , 

3~W  --  <*  +  S)F(x)- 

Combining  (4)  and  (5)  we  have  as  our  final  equation 
2 

~~ ^2^  ~  3ar(or  .+  S)F(x)  =  ^“(4«n2a  '  T(x)  4) 
dx  1  ’ 


(2.8) 


(2.9) 


We  wish  to  solve  this  equation  subject  to  the  boundary 
condition 


I,+(o)  =  p  I*  (o)  (2.10) 

where  p  is  the  mean  reflectance  for  diffuse  radiation  going  from 
the  material  into  the  vacuum. 

Now  the  equation  is  second  order  in  x  but  we  have  only  one 
boundary  condition.  However,  if 
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3o(o  +  s)  »  1, 


(2.11) 


we  expect  F{x)  does  not  change  rapidly  over  the  interior  of  the 
material  and  the  equation  is  of  boundary  layer  type.  In  any 
region  where 


d2F 

2 

« 

3o(o  +  S)F(x) 

dx 

* 

(2.12) 


then  we  have 


.  ~  4n2g 1  dT4 (x)  _  .  . 

F(x)  "  3(o  +  S)  dx  “  Fint  X  * 


(2.13) 


Near  x  =  0  this  interior  solution  is  modified  by  the  presence 
of  the  boundary.  Write 


F(x)  =  Fint(x)  +  G(x)  . 

Then  G(x)  satisfies 

,2  d2F.  .(x) 

—  ffl-  -  3 o(o  +  S)G(x)  =  -  - - 


dx 


dx 


(2.14) 


(2.15) 


But 


2 

d  F. 
int 

dx2 


will  be  small  so  we  simply  take 


2 

d-  -  3P  (ff  +  S)G(x) 

dx 


=  0. 


(2.16) 


Now  we  demand  that  G(x) 


0  rapidly  as 


0  /  J-...  ) 

LV  3o(o  +  S)  -1 


(2.17) 


So  G(x)  is  proportional  to 

exp.  [-  / 3o (o  +  s)  x] . 
We  discard  the  solution 

exp.  [+  / 3o (o  +  s)  x] 


5- 
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2.2  Spectral  Emission  of  Radiant  Energy 

In  the  preceding  section,  the  intensity  of  emitted 
radiation  I'“(o)  was  derived  as  a  function  of  the  total  in¬ 
tensity  of  forward  radiation  (i.e.,  the  Stefan-Boltzmann  Law 
was  used  to  derive  the  energy) .  This  derivation  assumed  that 
the  absorption  coefficient  a,  the  scattering  coefficient  S, 
and  the  index  of  refraction,  n,  were  constants.  This  is  not 
true,  however,  because  all  three  {a,  S,  and  n)  are  dependent 
upon  wavelength  (X) .  Appendix  I  characterizes  the  materials 
used  during  the  past  year;  the  data  tabulated  there  clearly 
indicates  that  the  three  properties  do  change  with  wavelength. 
In  addition,  the  absorption  coefficient,  or,  is  a  function  of 
temperature  as  well  as  wavelength. 

Because  of  this  wavelength  dependence  of  the  material 
properties,  either  Planck's  Law  or  Wien’s  Law  must  be  used  for 
the  calculation  of  the  intensity  of  spectral  emission.  Experi¬ 
mentally,  we  measure  the  energy  emitted  as  a  function  of  wave¬ 
length  rather  than  the  integrated  energy  over  all  wavelengths. 
Thus,  we  want  to  know  the  non- iso thermal  radiation  intensity 
as  a  function  of  wavelength  as  well  as  temperature  and 
temperature  gradient. 

The  exact  solution, when  the  temperature  gradient  is  known 
and  the  surface  temperature  must  be  determined, has  been  derived 
in  Appendix  II.  This  involves  the  solution  of  a  quartic  poly¬ 
nomial  equation  for  each  discrete  value  of  &  and  S.  To  avoid 
this  complication,  Wien's  Law  for  spectral  radiation, 

\  T  =  ^  "  C1  ^  5  exP(-c2/^T)  (2  *25) 

can  be  substituted  for  the  Stefan-Boltzmann  Law, 

I  =  n2g'T4  (2.26) 

Wien's  Law,  rather  than  Planck's  Law, 

c,X"5 

\,T  exp(+c2/XT)  -  1  (2.27) 

has  been  used  solely  for  the  reason  that  computations  are  more 
easily  carried  out.  The  error  introduced  by  using  Wien's  Law 
is  large  at  high  wavelengths  (25%  at  X  -  10|i  and  1000°  K)  but 
at  the  wavelengths  less  than  X  =  6p  where  scattering  is 
important  the  error  decreases.  For  example,  at  X  -  In  and 
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1000°K  it  is  only  0.000156,  and  at  X  ~  5d,  the  error  is  0.1%. 


The  substitution  of  Wien's  Law  directly  into  the  derived 
expression  of  the  previous  section  is  not  strictly  justified, 
in  the  sense  that  the  integration  of  Wien's  Law  over  all  wave¬ 
lengths  does  not  yield  the  Stefan-Boltzmann  energy  function. 

It  was  deemed  satisfactory  however,  because  the  region  of 
interest  for  non-iso thermal  and  isothermal  emissivity  as  a 
volume  process  in  ceramic  oxides  is  below  six  microns,  and  in 
this  region  the  errors  introduced  are  not  large.  If  we  remain 
within  the  region  l-5p,  the  error  is  less  than  one  percent. 


Substituting  Wien's  Law,  the  equation  (2.24) 

IMP)  : 


2 

n  a 


(^4cT"^)  +  s)  +  1//2  C1  +  P) 


["Tiff'd  +  - - -  +  - 1 -  i  y  24^ 

L  S3a{a  +"sj  dx  3flr(®  +  S)  ^2  J 


2  4 
T 

dx* 


becomes  ¥(o) 


n2CxX  5 

^  ^  S~3a{a  +  s)  +  1/2  (1  +  p) 


exP(-  7^)  + 


*-T  / 3<x  (&  +  s) 

C, 


d  6XP("  XT} 
dx 


+ 


2  ~2 

i _ d  sxp(-  s’ 


3a,(e>'  +  S) 


dx 


] 


(2.28) 
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or  f(o) 


r,2Clr5 


(X ,  T)  ,1  -  a,  _ _ 

(  4 a  )  ''3® (o'  +  s)  +  1/2(1  +  p) 


[-<-  S> +  ,i=  ft  -  <-  S’  £} + 


i 

3<*(ff  +  S) 


This  expression  includes  the  temperature  gradient  (™ )  , 

and  variations  of  the  absorption  coefficient  (©) ,  scattering 
coefficient  (S) ,  and  reflectance  (p)  with  wavelength.  and 
C2  are  the  first  and  second  radiation  constants,  X  the  wave¬ 
length,  and  T  the  absolute  temperature. 


Results  for  one  material  showing  the  ratio  of  non- iso¬ 
thermal  to  isothermal  energy  emitted  at  different  wavelengths 
are  plotted  in  Fig.  2.2.  This  ratio  is  a  measure  of  the 
expression 


,  2,  .  1  f  2  ,  2^  dTl 

exp(~  +  mx*'-*  W  ^  ^ 


(2.29) 


compared  to  exp(-  — ) .  Table  I  is  a  tabulation  of  the  specific 
values  used  in  computing  the  theoretical  curve. 


The  effect  of  increasing  the  magnitude  of  the  temperature 
gradient  is  to  increase  the  amplitude  of  the  curve.  The 
temperature  gradient  does  not  affect  the  wavelength  dependency 
of  the  energy  ratio. 

A  brief  description  of  the  three  black-body  radiation 
laws  is  given  in  Appendix  IV. 
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Figure  2.2  Ratio  of  Non-Isothermal  to  Isothermal 
Emissivity  for  AD-85  Calculated  from 
Equation  (2.28b). 
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TABLE  2-1 


Values  Used  to  Calculate  Non-I  sot  hernial  Energy  Ratio 
for  AD-85  Al20^  at  1100°C  Surface  Temperature 

\(p)  n,  o  or  (cm  3)  S(cm  1) 

K 


1 

1.75 

0.07 

0.513 

260 

0.051 

7.63 

X 

10-3 

2 

1.74 

0.07 

0.117 

260 

0.105 

3.62 

X 

io-3 

3 

1.72 

0.07 

0.256 

265 

0.070 

2.54 

X 

io~3 

4 

1.65 

0.07 

0.610 

263 

0.046 

1.91 

X 

io“3 

5 

1.60 

0.07 

5.420 

256 

0.015 

1.53 

X 

io-3 
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2.3  Calculation  of  Non-isothermal  Spectral  Eteissivities 

For  the  samples  of  alumina  and  strontium  titanate 
used  for  experimental  measurements,  non-isothermal  spectral 
emissivities  were  calculated  from  Equation  2.28b.  For  this 
calculation  values  of  the  absorption  coefficient,  scattering 
coefficient  and  reflectance  were  determined  as  described  in 
Reference  2.  Values  of  surface  temperature  and  temperature 
gradient  were  taken  to  correspond  with  the  experimental  con¬ 
ditions  actually  used.  Results  of  these  calculations  are  shown 
directly  compared  with  experimental  measurements  in  Figures 
3.5-3.24. 


These  results  for  the  calculated  values  of  the  ratio  of 
non-isothermal  to  isothermal  emissivity  indicate  that  above 
6  microns,  the  non-isothermal  effect  is  nil.  At  lower  wave¬ 
lengths,  the  effect  of  a  temperature  gradient  is  to  increase 
the  emissivity  whenever  there  is  a  positive  gradient,  and  to 
decrease  the  emissivity  whenever  there  is  a  negative  tempera¬ 
ture  gradient.  The  temperature  of  the  surface  itself  has  less 
effect  on  the  non-isothermal  emissivity  than  the  gradient. 

That  is,  the  term  containing  the  gradient  in  equation  (2.29) 
is  of  greatest  significance. 
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3.0  NON- ISOTHERMAL  EMISSIVITY  MEASUREMENTS 

3.1  Apparatus  for  Measuring  Non-Iso thermal  Emissivity 

The  purpose  of  the  apparatus  is  to  produce  a  con¬ 
trolled  temperature  gradient  across  a  disc-shaped  ceramic 
sample  whose  emissive  properties  are  to  be  measured.  The 
gradient  is  produced  only  in  the  radial  direction,  there¬ 
by  reducing  che  analysis  to  that  of  a  one  dimensional  heat 
transfer  problem.  A  rotating  sample  arrangement  is  used  to 
maintain  symmetric  temperature  conditions. 

The  apparatus  operates  inside  the  same  furnace  which  has 
been  used  for  isothermal  emissivity  measurement;  the  latter 
has  been  described  previously .  The  furnace  beating  elements 
are  used  to  adjust  the  surface  temperature.  An  internal 
resistance  element  is  used  to  produce  a  positive  temperature 
gradient  across  the  sample;  forced  air  cooling  through  a  hollow 
shaft  is  used  to  produce  a  negative  temperature  gradient  across 
the  sample.  The  complete  apparatus  is  shown  in  Figure  3.1. 

The  samples  and  the  internal  heating  element  are  supported 
on  a  one-half  inch  alumina  shaft.  The  alumina  shaft  is  sup¬ 
ported  and  rotated  by  stainless  steel  shafts  with  conical 
points  which  mate  with  conical  holes  at  each  end  of  the  alumina 
shaft.  The  two  stainless  steel  shafts  are  rotated  by  an 
electric  motor  and  synchronized  by  a  pair  of  seisyns.  This 
arrangement  assures  minimum  stress  and  vibration  on  the  shaft 
and  sample  assembly  while  giving  a  uniform  rate  of  rotation. 

In  order  to  minimize  the  effects  of  thermal  variations, 
the  apparatus  was  designed  to  produce  a  radial  temperature 
gradient  in  the  sample.  ine  gradient  is  obtained  in  a  sample 
which  has  a  width  of  one-half  inch  and  which  is  placed  between 
ceramic  discs  of  the  same  dimensions  and  made  of  a  material 
with  similar  thermal  properties.  The  discs  are  held  together 
by  an  axial  spring  load;  the  resulting  thin  air  gap  between 
adjacent  discs  reduces  axial  heat  transfer  and  resultant 
temperature  gradients  along  the  samples  and  support  plates. 

The  interior  heating  element  is  wound  for  a  length  of 
nine  inches  under  the  sample  and  support  plates,  with  the 
samples  located  at  the  center  of  the  winding.  This  arrange¬ 
ment  gives  a  uniform  axial  heat  zone  in  the  region  of  the 
sample  and  support  plate  assembly;  no  axial  gradient  is  pos¬ 
sible  since  there  is  no  variation  in  thermal  properties  in 
the  axial  direction. 
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Fumoce  shell  (water  cooled ) 


3.1.1  Internal  Heating  Element 

The  internal  heating  element,  producing  the  non- 
isothermal  condition  consists  of  a  Pt-40%  Rh  alloy  resistance 
element  rotating  with  the  samples.  This  wire  element  is 
spirally  wound  around  a  one-inch  O.D.  alumina  core.  The  core 
is  grooved  for  a  ten  turns  per  inch  winding  and  flame  sprayed 
with  alumina  to  cover  the  resistance  wire.  The  outer  surface 
is  ground  true  to  facilitate  assembly  of  the  sample  and  support 
plates  on  the  heating  element.  The  power  for  this  heater  is 
supplied  to  the  rotating  element  through  two  slip-ring 
assemblies  located  on  the  drive  shafts  at  either  end  of  the 
furnace . 


3.1.2  Internal  Cooling  Apparatus 

The  internal  cooling  of  the  emissivity  samples  is 
accomplished  by  forcing  compressed  air  through  the  hollow 
alumina  support  shaft  and  stainless  steel  shafts.  The  alumina 
shaft  is  sealed  to  the  stainless  steel  shafts  by  means  of  o- 
rings  and  set  screws  (Fig.  3.2) .  The  compressed  air  supply 
can  be  regulated  to  flow  air  through  the  shaft  at  from  0-15 
psi,  which  in  turn  controls  the  amount  of  heat  transferred 
away  from  the  specimen.  Gradients  which  have  been  established 
under  steady  state  conditions  for  alumina  samples  are  shown  in 
Fig.  3.3. 

Cracking  of  the  back-up  plates  when  subjected  to  an  in¬ 
crease  in  cooling  air  flow  while  at  temperature  was  a  major 
problem.  In  order  to  avoid  cracking,  it  was  necessary  to  flow 
the  cooling  air  at  the  chosen  rate  at  all  times.  This 
eliminated  the  complete  fracture  of  both  sample  and  back-up 
plates,  but  did  not  totally  eliminate  cracking. 

3.1.3  Measurement  of  Temperature  Gradient 

The  sample  temperature  and  temperature  distribution 
were  measured  by  a  series  of  thermocouples,  four  i.n  each  sample. 
The  thermocouple  junctions  were  arranged  in  a  radial  spiral 
covering  an  arc  of  90°  to  avoid  possible  thermal  interaction 
which  might  have  resulted  if  a  single  radial  line  were  used. 

Two  samples  were  normally  run  simultaneously,  requiring  one 
set  of  thermocouples  to  pass  completely  through  the  sample 
nearest  the  slip  ring  assembly. 

Thermocouple  elements  of  chromel-aluinel  sheathed  in  inconel 
tubing  or  platinum-rhodium  sheathed  in  platinum  tubing  were 
used.  For  temperatures  above  1000°C,  platinum-rhodium  elements 
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Alumina  Shaft. 


Sample  Temperature 


were  used  exclusively.  A  Leeds  and  Northrup  K-0  potentiometer 
was  used  to  measure  the  thermocouple  EMP.  A  low-noise  slip¬ 
ring  assembly  transferred  the  thermocouple  signal  from  the 
rotating  shaft  to  the  potentiometer. 

The  maximum  cuter  surface  temperature  achieved  under  non- 
isothermal  conditions  was  approximately  1100°C  (1373°K)  for 
both  negative  and  positive  temperature  gradients.  At  this 
outer  surface  temperature,  the  gradient  achieved  with  the 
inner  surface  at  a  higher  temperature  than  the  outer  surface 
(positive  gradient)  was  98  C°/cm.  The  negative  gradient 
achieved  at  this  surface  temperature  was  46C°/cm.  At  a  sur¬ 
face  temperature  around  905°C  (1178°K) ,  the  maximum  positive 
gradient  was  74  C°/cm  and  the  maximum  negative  gradient  was 
46  C°/cm.  The  maximum  positive  gradient  was  38  C°/cm  and  the 
maximum  negative  gradient  was  32  C°/cm  at  605°C  (878°K) . 

3.2  Emissivity  Measurement 

The  optical  system  used  for  non-isothermal  emissi¬ 
vity  measurement  is  shown  in  Pig.  3.4.  Input  to  the  spectro¬ 
meter  chopper  comes  (a)  from  the  sample  and  (b)  from  a  black- 
body  furnace  (described  in  reference  2)  maintained  at  a 
temperature  close  to  the  surface  temperature  of  the  sample. 

Because  the  blackbody  used  during  the  ncn-isothermal 
measurements  is  not  an  integral  part  of  the  emissivity  sample, 
temperature  control  and  corrections  for  temperature  differences 
become  critical.  A  difference  in  temperature  of  5°  between  two 
blackbodies  at.  1000°K  causes  difference  in  emitted  energy  of 
5%  at  a  wavelength  of  one  micron.  This  means  that  temperature 
must  be  accurately  controlled  as  discussed  in  reference  (2) . 
Otherwise,  emitted  energy  differences  due  to  inaccurate 
temperature  control  will  be  a  significant  contribution  to  the 
total  non-isothermal  radiant  energy  measured. 

In  the  measurements  reported  here,  temperature  differences 
between  the  blackbody  and  sample  surface  were  maintained  at 
less  than  2°c  (2%  error  at  X  =  1  n;  <  0.5;'  error  at  X  >  4  p)  . 

In  order  to  obtain  a  direct  measure  of  the  influence  of 
temperature  gradients  in  the  sample,  each  sample  was  measured 
under  (a)  isothermal,  (b)  positive  temperature  gradient,  and 
(c)  negative  temperature  gradient  at  the  same  surface  tempera¬ 
ture.  Results  were  then  analyzed  directly  as  the  ratio 
between  isothermal  and  non-isothermal  emissivity. 

In  addition  to  problems  of  temperature  measurement  there 
are  fluctuations  in  data  resulting  from  the  spectrometer 
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precision  of  0.5%  of  full  scale  reading.  Thus  as  discussed 
previously individual  measurements  are  estimated  as  having 
an  overall  precision  of  +  0.01  emissivity  units. 

3 . 3  Results 

Results  of  the  non-isothermai  emissivity  measure¬ 
ments  are  plotted  in  Figs.  3.5  -  3.24  as  the  ratio  of  non- 
isothermal  emissivity  (en)  to  isothermal  emissivity  (e^) 
versus  wavelength  (X)  .  On  these  figures,  calculated  values 
(section  2.3)  are  also  shown.  Measurements  were  run  at  sur¬ 
face  temperatures  of  approximately  600,  900,  and  1100 °C  for 
the  alumina  samples,  and  at  600  and  900°C  for  the  strontium 
titanate  and  Pyroceram  specimens.  The  maximum  gradient 
achieved  was  98C°/cm,  and  the  usual  gradient  30  to  50  C°/cm. 

3.3.1  Aluminum  Oxide 

Measured  and  calculated  values  of  the  emissivity 
ratio  for  the  four  alumina  specimens  are  shown  in  Fig.  3.5- 
3 . 16 . 


The  measured  curves  show  some  deviation  from  the  calcu¬ 
lated  values  over  the  region  1-5  microns,  but  the  measured 
ratio  is  within  the  estimated  accuracy  of  the  measurements. 

3.3.2  Strontium  Titanate 


Measured  and  calculated  results  for  strontium 
titanate  are  shown  in  Figs.  3.17  -  3.20.  The  experimental 
data  agrees  with  the  calculated  values  within  the  estimated 
accuracy  of  the  measurements. 


3.3.3  Pyroceram  9606  and  9608 

The  measured  values  for  Pyroceram  9606  and  9608 
*re  shown  in  Figs.  3.21  -  3.24.  These  curves  have  the  same 
general  shape  as  do  the  alumina  and  strontium  titanate  curves. 

3.4  Discussion 


The  spectral  region  of  most  interest  is  over  the 
wavelength  range  of  1  to  5  microns.  For  the  materials  studied, 
the  high  scattering  coefficients  lead  to  measured  isothermal 
emissivities  in  the  range  0.05  +  0.01  to  0.25  +  0.01.  At  these 
low  values  even  though  the  absolute  precision  is  good,  the 
relative  error  is  large  (20%  -  to  4%  of  the  measured  value) . 
Similar  results  are  found  for  the  non-isothermal  measurements. 
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Thus  when  we  compare  the  two  as  the  ratios  plotted  in  Pigs. 
3.5-3.24  the  relative  error  is  the  sum  of  the  two,  and  may 
be  as  large  as  +  40%  depending  on  the  absolute  level  of  the 
emissivities  involved.  The  most  common  absolute  values  are 
in  the  emissivity  range  0.15  -  0.20,  corresponding  to  a 
precision  of  +  15-10%  for  the  ratio. 

The  assumption  that  scattering  is  isotropic  in  the  for¬ 
ward  and  backward  directions  was  made  in  the  derivation  for 
the  intensity  of  emitted  radiation.  This  assumption  is  valid 
only  when  the  size  of  the  scattering  centers  is  small 
compared  to  the  wavelength  of  the  incoming  light,  or  for 
large  scattering  centers  compared  to  the  wavelength  of  the 
incoming  light  when  the  relative  index  of  refraction  (m)  of 
the  scattering  centers  is  zero  (Fig.  3.25).  For  large  pores 
in  alumina,  the  forward  scattering  is  approximately  three 
times  the  backward  scattering.  Therefore  the  measured  non- 
isothermal  to  isothermal  emissivity  ratio  is  expected  to  be 
somewhat  greater  than  the  calculated  value  for  a  positive 
temperature  gradient,  and  less  than  the  calculated  value  for 
a  negative  gradient.  For  strontium  titanate,  the  forward 
scattering  for  large  pores  is  about  one  and  one-half  times 
the  backward  scattering. 

In  the  wavelength  region  1-5  microns  where  emissivity 
is  a  volume  process,  there  is  agreement  between  the  measured 
and  calculated  values  within  the  limits  of  experimental 
precision.  In  all  cases  with  a  positive  temperature  gradient 
the  non-isothermal  emissivity  is  greater  than  the  isothermal 
value;  in  all  cases  with  a  negative  temperature  gradient  non- 
isothermal  emissivity  is  smaller  than  the  isothermal  value. 

The  differences  between  the  non-isothermal  and  isothermal 
values  are  clearly  significant,  having  ratios  as  large  as  1.35 
for  the  relatively  modest  temperature  gradients  imposed. 
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Figure  3.5  Emissivity  Ratio  vs.  Wavelength  for  AD-995 
at  approximately  600°C. 
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Figure  3.7  Emissivity  Ratio  vs.  Wavelength  for  AD-995 
at  approximately  1100°C. 
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Figure  3.8  Eraissivity  Ratio  vs.  Wavelength  for  AD-85 
at  approximately  600°C. 
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Emissivity  Ratio  vs.  Wavelength  for  AD-85 
at  approximately  900°C. 


-26- 


Ratio  (€n  /ex  ) 


X  {fM) 


Figure  3.10  Emissivity  Ratio  vs.  Wavelength  for  AD-85 
at  approximately  1100°C. 


-cl- 


Ratio  (6n/€x 


X  (p) 


Figure  3.11  Emissivity  Ratio  vs.  Wavelength  for  &£-l 
at  approximately  600 °C- 
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Figure  3.12  Emissivity  Ratio  vs.  Wavelength  for  AjJ-l 
at  approximately  900°C. 
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Figure  3.16  Emissivity  Ratio  vs.  Wavelength  for  hi  4 
at  approximately  1100°C- 
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Figure  3.17  Emissivity  Ratio  vs.  Wavelength  for  SrTiC^A 
at  approximately  600°C. 
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Figure  3.19  Emissivity  Ratio  vs.  Wavelength  for  SrTi03B 
at  approximately  600°C. 


-36- 


Figure  3.20  Emissivity  Ratio  vs.  Wavelength  for  SrTiC^B 
at  approximately  900 °C- 
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Figure  3.21  Emissivity  Ratio  vs.  Wavelength  for  Pyroceram 
9606  at  approximately  600°C. 


-38- 


€n/€ 


Rq\\0(€^/€1) 


1.5 


1*4|  O  Positive  Gr 


0  2  4  6  8 


X  {/J,) 


Figure  3.24  Emissivity  Ratio  vs.  Wavelet 
9608  at  approximately  900°C. 
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Figure  3-25  Relative  Fraction  of  Back  Scattering  for  Pores 
of  Large  Size  in  a  Solid  Matrix. 


4.0  TRANSMISSION  OF  RADIANT  ENERGY  THROUGH  A  MATERIAL  WITH 
ABSORPTION  AND  SCATTERING 


When  a  beam  of  parallel  light  is  incident  normally  upon 
a  thin  slab  containing  scattering  centers,  the  emerging  light 
may  consist  of  both  parallel  and  diffuse  light.  The  differ¬ 
ential  equations  describing  this  system  have  been  formulated 
and  solved  by  Ryde(3)  .  Some  of  Ryde's  equations  have  been 
reproduced  here  to  have  them  available  in  the  notation  used 
in  this  report . 


The  equations  for  the  forward  and  backward  diffuse  flux 


are: 


j  T 

—  =  f'l  -  (a  +  s)  I  +  sj 
dx  P 


(4.1) 


-  ~  =  s'l  -  (a  +  s)J  +  si 
dx  P 

and  for  the  forward  and  backward  parallel  flux: 
d! 

=s  -  (a  +  s'  +  f)I  =  -  g'l 
dx  p  p 


dJ 


dx 


5  -  -  ( 


a+s  +  f  )  J  =  -  g  J 
P  P 


(4.2) 


(4.3) 


(4.4) 


.(3) 


The  solutions  to  these  equations  obtained  by  Ryde  are: 
I  =  K  exp(-g'x)  (4.5) 


J  =  M  exp(g  x) 
P 


(4.6) 


I  =  A(i  -]3  )  exp(aQx)  +  B(1  +  $o)exp(-aQx)  - 


U  exp(-g'x) 


(4.7) 


where  U 


J  =  A(1  +  j3Q)  exp  (aQx)  +  B(1 
V  exp(-g'x) 

_  . +_s  f  g'l 

.2  2 

g  -  a 

o 


0  ) exp(-cQx) 


(4.8) 

(4.9) 
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V 


(4.10) 


sf*  +  s* (a  +  s  -  q* ) 

,  2  2 
g  -  ao 

and  K,  M,  A,  and  B  are  determined  by  the  boundary  conditions. 

If  the  amount  of  paraJlel  light  passing  through  the 
specimen  is  large  (i.e.,  very  thin  sample  and  very  small 
scattering  coefficient) ,  then  equations  4.1  and  4.2  are  not 
correct,  but  must  contain  a  term  representing  the  sc<  tturing 
from  the  parallel  flux  in  the  negative  x-direction.  This  flux 
in  the  (-x)  direction  i. v  caused  by  reflection  from  the  sample- 
air  interface  where  the  beam  is  emerging  from  the  sample. 

These  conditions  and  the  differential  equations  are  treated  in 
Appendix  III. 

For  the  case  where  the  sample  is  sufficiently  thick  and 
the  scattering  coefficient  large,  the  intensity  of  the  trans¬ 
mitted  diffuse  and.  parallel  light  beams  may  be  solved  as 
follows: 

1)  Consider  the  case  wherein  the  index  of  refraction 
of  the  medium  containing  the  scattering  centers  is 
the  same  as  the  index  outside  the  medium.  Then 
there  are  no  reflections,  and  the  boundary 
conditions  for  diffuse  light  for  a  beam  of  parallel 
light  of  unit  intensity  incident  normally  on  the 
medium  are 

1)  at  x  =  o,  I  =  o,  J  =  R1  (4.11) 

2)  at  x  =  d,  I  =  T1  -  I  ,  J  =  o  (4.12) 

where  T1  is  the  total  light  transmitted  both  para¬ 
llel  and  diffuse  and  R1  is  the  amount  reflected. 

Then  ?.  system  of  four  equations  in  terms  of  A,  B, 

T1 ,  and  R1  may  be  solved  for  A,  B,  T1 ,  and  R1 . 

These  equations  are: 

A  (1  -  3  )  +  B  ( 1  +  3  )  =  U  (4.13) 

A ( 1  +  3  )  +  B ( 1  -  3  )-  R1  =  V  (4.14) 
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(4.15) 


A(1  -  0  )exp(o  d)  +  B(1  +  0  )exp(-a  d)  -  T*  * 
o  o  o  o 

(U  -  l)exp(-g'd) 

A(1  +  0Q)exp(aod)  +  B(1  -  0o)exp(-aQd)  = 

Vexp(-g'd)  (4.16) 

Therefore: 

t  =  2(1  +  0^)  sinh  a  d  +  40  cosh  a  d  (4.17) 

o  o  o  o 

A  =  ~  {u ( 1  -  0o)exp(-aQd)  -  V(1  +  0q) exp(-g'd) }  (4.18) 

B  =  ~  {V(l  -  0q) exp(-g 1 d)  -  U(1  +  0o)exp(aod)}  (4.19) 

It  is  not  necessary  to  solve  for  t‘  and  R1 ,  since 
equation  (4.6)  does  not  contain  either  R1  or  T1 . 
Substituting  A  and  B  into  equation  (4.7)-  (4.7) 

Ix=d  C(1  -  @Q)  Ud  -  30)exp(-tJod)exp((Jod)  - 

(1  -  0q)  V{1  +  0q) exp(-g'd) exp(aQd)  + 

(1  +  0  )  V (1  -  0  ) exp(-g'd)  exp (-ad)  - 
00  o 

(1  +  0q)  U(1  +  0q) exp(aQd)exp(-aod) }  - 

Uexp(-g'd)  (4.20) 

I  .  =  T  (40  U  +  V(1  -  0^)exp(-g'd)  sinh  a  d}  - 
x=d  t  o  o  o 

Uexp(-g'd)  (4.21) 

and  for  parallel  light 

I  =  exp(-g'd)  (4.22) 

Px=d 
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Next,  consider  the  case  where  there  are  reflections 
from  the  surfaces  of  the  specimen.  For  the  parallel 
segment,  the  boundary  conditions  are  (Fig.  4.1): 


1)  x  —  d ,  J  *  I  p  ,  »  0 

P  Jr 


(4.23) 


2)  x  *  o,  I  ••  (1  -  p  ) 
p  '  Kno 


(4.24) 


c_  =  (1  -  pnQ) exp(-g'd) 
px=d  ^ 


(4.25) 


For  the  diffuse  portion: 

I  -  ■  (1  -  p  )  [7  {40  U  +  V(1  -  02)  exp(-g'd)  sinh  and) 
x=d  ^no  t  o  o  ^  3  0 


Uexp(-g*d) 


(4.26) 


The  intensity  of  the  light  emerging  from  the  surface 


1)  For  parallel  light  I1 


*  (1  -  D  .)I 
ni  p 


(1  -  p  .)  (1  -  p  )exp(-g  d) 
pni  yno  *  ^ 


(4.27) 


or  I  =  (1  -  p  .;  exp ( -g  d) 
p  ,  'nx 

x=d 


(4.28) 


2)  For  diffuse  light  =  (1  - 

or  I*  =  (1  -  p .)  (1  -  p  )  [7  {40  U  + 
x=a  Kn o'  '‘t  o 


(4.29) 


V ( 1  -  0^)exp(-g  d) sinh  aQd}  - 


Uexo 


(-g‘ a) ] 


(4.30) 


Since  the  intensity  of  the  incident  light  was  taken 
as  unity,  I*  and  I*  are  the  transmissivities, 

r  j 

x=d 

Tp  and  of  the  material  for  parallel  and  diffuse 
light . 
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Absorbing  and 
scattering  material 


x  =  o  x  =  d 


Fig.  4.1  Heat  flew  by  combined  radiation  and  conduction 
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The  materials  investigated  under  this  cor t -act  all  have 
large  scattering  coefficients  (see  Appendix  A) ,  For  a 
specimen  thickness  of  0.006"  (0.015  cm)  and  scattering 
coefficients  s'  and  f'  of  90  and  260  cm”!,  with  a  =0.8  and 
pni  “0.07  (corresponding  to  AD-85  AI2O3  @  X  *  lp  and  T  = 
500°C) ,  a  calculation  of  I*  reveals  that: 

px-d 

X  =  (0.93)exp(-5.26)  =  0.005  (4.31) 

px>d 

X*  «  (0.93) 2exp(-5.26)  =  0.004  (4.32) 

Px=d 

or  a  transmission  for  parallel  light  of  less  than  one  percent. 
From  equation  4.23, 

J  =  (0.07)  (0.005)  =  0.0004  (4.33) 

Px=d 

which  justifies  neglecting  the  scattering  from  the  backward 
parallel  flux  in  calculating  the  diffuse  flux.  Equations  (4. 
28)  and  (4.30)  have  been  used  to  calculate  the  transmission 
of  the  thin  sections  described  in  the  next  section. 


5.0  THIN  SAMPLE  TRANSMISSION  MEASUREMENTS 


Thin  sample  transmission  measurements,  when  made  on 
single  crystal  specimens,  allow  evaluation  of  the  absorption 
coefficient  (<*)  of  the  single  crystal.  The  equation 

=  (1  -  P_.)2exp(-g,d)  (5.1) 

y  px=d 

allows  evaluation  of  (g‘),  the  extinction  coefficient  for  poly¬ 
crystalline  thin  sections.  This  extinction  coefficient  is  the 
sum  of  twice  the  single  crystal  absorption  coefficient  plus  the 
forward  and  backward  scattering  coefficients  for  parallel  light. 

5.1  Apparatus 

It  was  originally  planned  that  the  single  crystal 
transmission  furnace  and  optics (1)  be  used  so  that  measurements 
could  be  taken  at  elevated  temperatures.  The  low  transmission, 
of  energy  through  the  sample  prevented  the  use  of  both  che 
furnace  and  the  special  optics.  Instead,  the  sample  was  in¬ 
serted  into  the  normal  optical  path  of  the  spectrophotometer 
(Fig.  5.1) .  Room  temperature  measurements  only  were  run 
because  there  was  no  furnace  which  was  small  enough  to  fit  into 
the  limited  space  available  on  the  spectrophotometer. 

5.2  Transmission  Measurements 

The  amount  of  parallel  energy  passed  through  the 
sample  was  very  small.  As  calculated  in  section  4,  the  amount 
of  parallel  light  (wavelength  l^1)  passed  by  a  .006"  sample  of 
AD-85  at  500°C  is  0.4%.  Thus,  measurements  were  being  taken 
at  the  extreme  limits  of  the  spectrophotometer  precision. 

In  order  to  obtain  measurements,  the  following  technique 
was  used: 

1.  A  curve  was  run  with  both  optical  paths  open 
(referred  to  as  "100%  line") . 

2.  A  curve  was  run  with  the  sample  path  completely 
blocked  ("zero  percent"  line) . 

3.  A  curve  was  run  with  the  sample  inserted  into  the 
sample  optical  path. 
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4. 


The  reference  beam  was  partially  blocked  so  that  the 
energy  ratio  was  smaller. 

5.  Th>i  sample  curve  was  re-run  at  this  amplifier 
setting . 

6.  The  curves  (1,  2,  3,  and  5)  were  normalized  to  give 
transmission  data  for  the  thin  sections. 

This  process  was  tried  using  the  single  crystal  furnace  and 
optical  system,  but  was  unsuccessful  because  of  energy  losses 
by  the  optics. 

For  diffuse  tran.f  nission  measurements  it  is  desirable  to 
measure  the  transmitt  1  radiation  over  a  solid  angle  of  180°; 
for  parallel  transmission  measurements  the  ideal  arrangement 
would  be  a  collimated  measurement  system  assuring  that  only 
the  parallel  beam  reaches  the  instrument  collector.  In  the 
present  measurements  the  optical  system  collects  transmitted 
radiation  on  a  solid  angle  of  15°  fixed  by  the  minimum  energy 
input  necessary  for  barely  acceptable  precision.  Thus,  the 
intensities  measured  are  that  of  the  parallel  radiation  plus 
the  diffusa  radiation  falling  within  a  15°  solid  angle  to  the 
parallel  axis. 

5.3  Results 

Transmission  curves  were  run  on  the  alumina  samples, 
two  strontium  titanate  samples,  and  two  Pyroceram  samples. 
These  results  (extinction  coefficient  vs.  wavelength)  are 
plotted  in  Figs.  5.2  to  5.8.  Two  samples  of  alumina  AD-995 
were  run,  one  sample  being  0.006"  thick  and  the  other  specimen 
0.004"  thick.  The  results  (Fig.  5.2)  indicate  that  experimen¬ 
tally  the  extinction  coefficient  is  not  the  same  for  both 
samples.  Physically,  the  extinction  coefficient  should  be 
independent  of  thickness.  The  difference  between  the 
extinction  coefficients  is  thus  an  indication  of  the  magnitude 

of  the  errors  involved  in  the  transmission  measurements. 
Deviations  between  the  two  curves  range  from  12%  to  84%. 

For  the  AD-995  specimens,  the  transmission  at  one  micron 
for  the  0.006"  specimen  is  0.0507,  or  5.07%,  and  for  the 
0.004"  specimen  it  is  0.0568,  or  5.68%.  Data  on  the  spectro¬ 
photometer  (4)  give  an  accuracy  of  0.5%  3nd  a  reproducibility 
of  0.25%  for  the  ordinate  scale  of  the  instrument.  This  would 
mean  an  estimated  maximum  error  for  the  instrument  of  0.50% 
of  full  scale  when  comparing  two  specimens  of  the  same 
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Figure  5.3  Extinction  Coefficient  vs.  Wavelength  for  AD-85. 
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Figure  5.7  Extinction  Coefficient  vs.  Wavelength  for 
Pyrocerain  ?6  >b. 
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material  under  the  same  conditions.  For  energy  transmissions 
of  less  than  15%  (full  scale)  ,  this  would  be  minimum  error 
of  3%,  and  could  be  50%  or  greater  whenever  the  transmission 
is  less  than  1%. 

The  properties  of  the  materials  tested  were  such  that  we 
were  talcing  measurements  in  the  region  of  15%  or  less  for  the 
total  parallel  transmitted  energy.  Therefore,  it  does  not 
seem  unreasonable  to  assume  an  error  of  50%  for  the  measured 
coefficient  of  extinction  of  the  polycrystalline  specimens. 

5.4  Comparison  with  Calculated  Results 

Values  of  the  diffuse  transmissivity  and  parallel 
rransmissivity  were  calculated  from  equations  (4.28)  and  (4.30) , 
using  the  IBM  7090  computer  at  The  Aeronautical  Systems 
Division,  Wright-Patterson  Air  Force  Base,  Ohio. 

The  extinction  coefficient  g*  was  calculated  from  material 
properties.  By  definition  (equation  4.3  and  4.4) 

g'  =  a  +  s  +  f  (5.2) 

where  a,  s',  and  f*  are  as  previously  defined  (see  equation 
5.1) .  The  normal  scattering  coefficient,  S,  bears  the 
following  relationships  to  the  forward  and  backward  scattering 
coefficients  for  parallel  light,  f  and  s‘  respectively; 

S  =  f'  +  s'  (5.3) 

(2) 

And  from  Folweiler  and  Mallio  ,  for  alumina 

f*  ^  3  s'  (5.4) 

The  scattering  coefficient,  S,  is  related  to  the  number 
of  pores  per  unit  volume,  N,  and  the  pore  radius,  r,  and  the 
scattering  factor,  K,  by  the  relation 

S  =  KN(Trr^)  (5.5) 

If  the  volume  fraction  porosity  is  P,  then 


The  scattering  factor  K  varies  between  zero  and  four  and 
depends  on  the  pore  radius,  r,  the  wavelength  of  radiation. 
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9.0 


Values  of  the  Scattering  Factor  K  vs .  p  for 
Various  Values  of  the  Relative  Index  m,  when 


X,  and  the  relative  index  of  the  scattering  center  and  medium, 

m. 


The  scattering  factor,  K,  has  been  summarized  as  a 
function  of  the  dimensionless  parameter, 

P  =  —r1 1  (s.7) 

by  Van  de  Hulst^l  as  shown  in  Fig.  5.9,  for  values  of  m  <  1. 


Using  equation  (5.6)  and  the  relation  that  a,  the  absorp¬ 
tion  coefficient  appropriate  to  the  Ryde  equation,  is  equal  to 
twice  <*,  the  single  crystal  absorption  coefficient,  it  is 
possible  to  compute  the  extinction  coefficient. 


Calculation  of  values  of  the  extinction  coefficient  were 
included  in  the  computer  program  to  evaluate  the  expression 
for  the  diffuse  transmission  (equation  4.29) .  These  calcu¬ 
lations  indicate  that  the  extinction  coefficients  are  greater 
than  the  values  measured  with  the  spectrophotometer.  Table 
5-1  is  a  tabulation  of  calculated  and  measured  extinction 
coefficients  for  alumina,  strontium  titanate,  and  Pyroceram 
samples . 


For  most  of  the  sample  conditions  investigated,  it  was 
sufficiently  accurate  to  set  the  extinction  coefficient,  g' , 
equal  to  (s‘  +  f'),  or  S,  because  the  absorption  coefficient 
a,  was  very  small  compared  to  S,  especially  in  the  region 
1-4  microns.  It  may  be  concluded,  therefore,  that  the  anni¬ 
hilation  of  a  beam  of  parallel  light  by  a  body  having  an 
absorption  constant  a  <  <  s,  the  scattering  coefficient,  is 
caused  by  scattering. 

If  the  extinction  coefficient  g'  were  not  great  enough 
or  the  sample  was  not  thick  enough  to  neglect  the  parallel 
beam,  difficulties  about  the  uniformity  of  the  diffuse  beam 
would  arise.  That  is,  the  scattered  light  would  not  be 
homogeneously  dispersed  in  either  the  forward  or  backward 
direction. 


The  measured  transmissivity,  and  the  calculated  parallel 
and  diffuse  transmissivities  for  the  alumina  and  strontium 
titanate  samples  are  shown  in  Figs.  5.10  -  5.15.  As  can  be 
seen,  the  measured  transmissivities  are  greater  than  the 
parallel  transmissivity  (equation  5.1),  but  less  than  the 
diffuse  transmissivity  (equation  4.29).  The  measured  trans¬ 
missivity  approaches  the  total  transmissivity  (t^  +  t  )  as 
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ssivity 


Figure  5.10  Measured  and  Calculated  Transmissivity  vs. 
Wavelength  for  AD-995. 
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Figure  5.11  Measured  and  Calculated  Transmissivity  vs. 
Wavelength  for  AD-85. 
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calculated  diffuse  transmissivity 


0.12 


Wavelength  (microns) 


Figure  5.12  Measured  and  Calculated  Transmissivity 
Wavelength  for  Ai-1. 


Trans 


Figure  5.13  Measured  and  Calculated  Transmissivity  vs. 
Wavelength  for  Ai-4. 
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Figure  5.14  Measured  and  Calculated  Transmissivity  vs. 
Wavelength  for  SrTiO^A. 
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Transmissivity 


Wavelength  (microns) 


Figure  5.15:  Measured  Transmissivity  vs.  Wavelength 
for  Pyroceram  9605 . 
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soivity 


Wavelength  (microns) 


Figure  5.16:  Measured  Transmissivity  vs.  wavelength 
for  Pyrocerara  S608. 
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the  vav^ler  **_h  Increases.  This  means  that  the  angular  dis¬ 
tribution  of  transmitted  light  is  not  ideally  diffuse  (i-e., 
does  not  follow  a  cosine  law) .  but  has  an  angular  distribu¬ 
tion  in  the  forward  direction  such  that  more  is  transmitted 
in  the  area  subtended  by  the  optical  system  than  in  the  high 
angle  region.  The  infrared  scattering  work  of  Love  and 
Wheasler(6) .  using  an  aerosol  of  six  micron  alumina-particles 
as  the  scattering  medium  show  that  the  amount  of  the 
scattered  light  in  the  angular  range  0-10  degrees  is  60  to 
85%  of  the  total  amount  scattered  in  the  forward  direction, 
depending  on  the  wavelength  of  the  incoming  light.  The 
theoretical  work  of  Mie  (see  Ref.  (5) ,  Appendix  E)  indicated 
that  as  a  function  of  wavelength,  tha  intensity  at  low  angles 
(0-20°)  decreases  with  the  wavelength  when  the  scattering 
media  contains  absorbing  particles.  The  opposite  is  true  for 
our  system  of  a  scattering  media  with  non-absorbing  particles - 
This  difference  in  the  angular  intensity  distribution  with 
wavelength  is  slight,  however. 
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TABLE  5-1:  Calculated  and  Measured  Extinction  Coefficient*  (cm 
All  samples  6  mils  thick  except  where  noted. 
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♦Polished  specimens  have  metallographic  polish  cn  surfaces. 


6.0  COMPA'ISON  BETWEEN  THE  EXACT  AND  HA MAKER  SOLUTIONS  FOR 
HEAT  fRANSMISSION  BY  COMBINED  CONDUCTIVE  AND  RADIATIVE 
TRANSFER 


The  exact  equations  for  describing  the  transmission  cf 
heat  by  combined  conductive  and  radiative  transfer  in  a  homo¬ 
geneous  medium  are  (for  the  one -dimensional  case  and  the 
steady  state) : 

2 

k  -  2  a  n2  c*  T*  +  a (I  +  J)  =  0  (6.1) 

ax 


dl  2.4 

—  =  an  c  T  -(a  +  s)I  +  sJ 


(6.2) 


dj 

dx 


(a  +  s)  J  -  si 


2 

a  n 


(6.3) 


where  I  and  J  are  the  radiant  fluxes  in  the  +x  and  -x  direction, 
T  is  the  absolute  temperature,  c*  is  the  Boltzmann  constant, 
n  is  the  refractive  index,  k  is  the  thermal  conductivity,  a 
is  the  absorption  coefficient  for  diffuse  radiation,  and  s 
is  the  back-scattering  coefficient  for  diffuse  radiation. 


It  is  convenient  to  reduce  the  set  of  Eq.  (6.1),  (6.2), 

and  (6.3),  to  a  single  second-order  equation.  Eq.  (6.2)  and 
Eq.  (6-3)  are  added  and  the  sum  integrated  once  to  give: 

I  +  J  =  -(a  +  2s)  (H  x  +  kT)  +  A  (6.4) 

in  which  A  is  a  constant  of  integration  and  H  is  the  constant 
heat  flux  defined  by: 

dT 

H  =  I  -  J  -  k  —  (6.5) 

Substituting  in  Eq.  6.1  we  obtain: 

2 

k  r  -  2  a  n  a  T4  -  a  (a  +  2s)  (Hx  +  kT)  +  a  A  =  0  (6.6) 
dx^ 


For  simplicity  we  consider  the  one -dimensional  flow  of 
heat  by  conduction  and  radiation  inside  a  medium  as  shown  in 
Fig.  4.1.  Once  the  material  constants,  k,  a,  s,  and  n^  are 
specified,  there  remain  four  independent  initial  conditions 
to  be  fixed;  i.e.,  TQ,  (dT/dx)Q,  IQ,  and  JQ.  For  various 
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values  of  these  initial  conditons  we  have  constructed  a  family 
of  solutions  by  numerical  integration  of  Eq.  (6.6) - 

The  calculation  consists  of  (1)  selection  of  material  con¬ 
stants  and  initial  conditions;  (2)  evaluation  of  A  and  H  from 
Eqs.  (6.4)  arid  (5.S);  (3)  calculation  of  T(x)  by  numerical  in¬ 
tegration  of  Eq-  (6.6);  (4)  calculation  of  k(dT/dx)  ;  a::d  (5) 
calculation  of  I(x)  and  J(x)  from  the  relations: 

I  *  [-(a  +  2s)  (Bx  +  kT)  +  a  +  H  +  M^)]  (6.7) 

4  ax 

i  dT 

J  =  T  [-U  +  2s)  (Kx  +  kT)  +  C  -H  -  k (— — )  ]  (6.8) 

4  GX 

Bamaker's  approximation  consists  of  replacing  the  first 
term  in  Eq.  (6.1)  k(d2T/dx2)  by  k/Tg  d2(T4)/dx2  in  which  TQ  is 
some  average  temperature  in  the  material.  This  replacement 
makes  the  set  of  Eqs.  (6.1),  (6.2),  and  (6.3)  linear  and  homo¬ 
geneous  in  the  quantities  I,  J,  and  T^  so  that  a  solution  in 
closed  form  can  be  obtained  readily. 

The  Hamaker  solutions  for  the  radiant  fluxes  are  given  in 
terms  of  the  lumped  material  constants: 

b  =  4n2  o'  T3  (6.9) 


a  =  V  ^  +  a(a  +  2s)  =  ac  /(l  +  *>  (6.10) 

8  =  a/  (a  +  2s)  (6.11) 

k  =  2b/k  (a  +  2s)  =  2bSAa  (6.12) 

Then,  the  radiant  energy  fluxes  and  temperature  are  given  by: 
I  =  A ( 1  -  8)  e°X  +  B ( 1  +  8)  e-CTX  +  C(ox  -  B)  +  D  (6.13) 

J  =  A(1  +  B)  eCTX  +  B(1  -  8)  e~UX  +  C(ax  +  8)  +  D  (6.14) 

ox  —ox 

E  =  -AKe  -  B/ce  +  Ccx  +  D  (6.15) 


(6.16) 
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where 


J  0(1  +  3)  -  E  (23o>  -  (|f  )  20 
o  o  dx 

_ o 

4 0  a(* c  +  1) 

-J  0(1-3)  -  E  ( 23o)  +  (f§-)23 
o  o  ax 

_ o 

40  a(#c  +  1) 

nP 

J  2/c a  +  (5^  )  20 
o  dx 

_ o _ 

40  a(/c  +  1) 

j  2 3fCC  +  E  (40a) 
o  o 

43a  (>c  +  1) 


(6.17) 


(6.18) 


(6.19) 


(6.20) 


The  calculation  consists  of  (1)  selection  of  material 
constants  and  initial  conditions;  (2)  evaluation  of  b,  a,  0, 
and  ■<  from  Eq.  (6.9),  (6.10),  (6.11),  and  (6.12);  (3)  evalua¬ 
tion  of  A,  B,  C,  and  D  from  Eqs.  (6.17),  (6. 18),  (6.19),  and 
(6.20);  and  (4)  calculation  of  I,  J,  E,  and  T  for  various 
values  of  x  from  Eqs.  (6.13),  (6.14),  (6.15),  and  (6.16). 


A  comparison  has  been  made  between  numerical  integration 
of  the  exact  equation  and  numerical  evaluation  of  the  Hamaker 
solution  for  the  system  described  for  a  range  of  initial  con¬ 
ditions  and  material  constants  corresponding  to  previously 
measured  values  of  a ,  s,  and  n.  Computations  were  done  on  an 
IBM  7090  computer  at  the  Aeronautical  Systems  Division,  Wright 
Patterson  Air  Force  Base. 


For  all  calculations  IQ  was  set  equal  to  zero;  physically 
this  corresponds  to  radiation  from  the  surface  into  a  semi- 
infinite  media  with  the  same  index  as  the  test  section  (no 
surface  reflection)  and  having  perfect  transmissivity  (no 
scattering  or  absorption) .  Values  of  JQ  ranged  from  0.26-125 
cal  cnT^  sec-^,  values  of  (dT/dx)  were  100°K  cm-^  and  values 
of  TQ  were  1000  and  2000 °K.  Material  constants  were  selected 
to  correspond  to  measured  values  of  a ,  s,  n  and  K  for  A^O^, 
MgO,  S.i02  and  SrTi03.  Tabulated  input  data  and  a  comparison 
of  results  between  the  exact  and  Hamaker  solutions  are  given 
in  Table  6.1, and  Figs.  6.1  and  6.2,  are  typical  results. 

The  depth  of  the  sample  contributing  to  the  emissivity 
characteristics  is  that  part  that  is  optically  "near"  the 
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Figure  6.1  Comparison  Between  Exact  and  Hamaker  Solution 
Using  Conditions  Approximating  Test  Material 
Properties,  k  =  0.020;  a  =  0.23;  S  =  260; 

?  dm 

n  =  3;  (~)  =  100;  T  =  1,000;  J  =  5;  I  « 
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surface.  Thus,  even  though  calculations  were  carried  out  to 
a  greater  sample  depth  in  most  cases,  we  have  compared  results 
at  x  =  (a  +  2s)  and  x  *=  2  (a  +  2s)  “1  in  Table  6.1.  As  shown 
there  for  the  entire  range  of  conditions  selected  agreement 
between  the  exact  and  the  Hamaker  equations  was  within  0.1%. 
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7.0  SUMMARY 


At  wavelengths  in  the  range  from  one  to  six  microns , oxide 
materials  are  at  least  partially  transparent.  As  a  result 
thermal  emissivity  is  a  volume  process  dependent  on  the 
sample's  radiant  energy  absorption  and  scattering  character¬ 
istics.  In  Parts  and  II (2)  ©f  the  present  study, the  in¬ 

fluence  of  these  material  and  specimen  properties  in  isothermal 
emissivity  was  evaluated  analytically  and  measured  experi¬ 
mentally.  Results  showed  that  the  behavior  of  these  materials 
could  be  analyzed  satisfactorily  by  means  of  the  Hamaker 
approximation,  and  that  calculations  based  on  optical  character¬ 
istics  and  microstructure  features  were  in  satisfactory  agree¬ 
ment  with  experimental  measurements. 

In  th°  present  report ,we  have  extended  these  analyses  to 
specifically  deal  with  the  question  of  radiant  energy  emission 
from  a  non-isothermal  body,  in  Section  2 . 0 , analytical  evalu¬ 
ations  of  the  problem  for  total  energy  emission  (assuming 
material  constants  independent  of  wavelength)  and  for  the  more 
important  problem  of  spectral  radiant  energy  emission  (which 
allows  one  to  account  for  the  variation  of  refractive  index, 
absorption  coefficient  and  scattering  coefficient)  are  reported. 
On  the  basis  of  this  development  it  is  possible  to  calculate, 
and  we  have  calculated,  the  influence  of  temperature  gradients 
on  radiant  energy  emission. 

In  Section  3.0, a  device  for  measurement  of  non-isothermal 
emissivity  has  been  described  and  experimental  results  are 
compared  with  our  calculations.  This  comparison  is  based  on 
individual  measurements  precise  to  0.01  emissivity  units.  At 
the  wavelengths  of  most  interest  we  are  concerned  with  low 
emissivity  values-in  the  range  0.05-0.25.  As  a  resuxt,  even 
with  this  small  absolute  error,  the  precision  of  the  comparison 
is  not  high.  However,  calculated  and  measured  values  agree 
within  our  rather  broad  error  limits.  The  methods  of  material 
characterization  and  emissivity  evaluation  presented  are  clearly 
applicable  to  non-isothermal  as  well  as  isothermal  conditions. 

As  predicted,  the  influence  of  temperature  gradients  on  radiant 
energy  emission  is  found  to  be  substantial  and  of  real 
importance  for  heat  transfer  calculations. 

Thus, it  is  confirmed  that  in  this  1-5  micron  range  raaiant 
energy  emissivity  can  in  no  way  be  regarded  as  a  material,  or 
even  as  a  sample,  characteristic.  It  is  influenced  by  the 
microstructure  (mostly  porosity) ;  but  it  is  also  influenced  by 
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the  temperature  gradients  in  the  solid  as  well  as  the  surface 
temperature-  To  a  certain  extent  one  can  academically  cir¬ 
cumvent  this  difficulty  by  always  including  the  word  “iso¬ 
thermal”  in  definitions,  but  this  merely  avoids  the  problem. 

In  the  present  work  we  have  unequivocably  shown  that  the  effect 
of  temperature  gradients  is  large,  and  we. have  tried  to  show- 
how  to  cope  with  it. 

Calculations  of  the  transmission  of  radiant  energy  are 
presented  in  Section  4.0  of  the  present  report-  Comparisons 
of  these  calculations  with  experimental  measurements  (Section 
5.0)  offer  still  greater  problems  than  for  emissivity  because 
of  the  low  values  of  transmission  we  are  dealing  with  and  con¬ 
sequent  low  precision.  Within  the  limits  of  experimental 
precision,  calculations  and  experiment  are  in  agreement.  How¬ 
ever  the  difficulty  of  measuring  aU  the  diffuse  transmission 
(through  a  solid  angle  of  180°}  and  the  extremely  low  fraction 
of  parallel  transmission  make  an  exact  analysis  of  the  parallel 
and  diffuse  components,  and  a  clear  separation  of  the  two, 
unattainable  with  the  apparatus  and  techniques  employed  here- 

It  was  hoped  that  evaluation  of  the  parallel  and  diffuse 
transmission  would  lead  to  a  basis  for  determining  the 
effective  absorption  coefficient  in  oxide  samples  as  a 
separable  variable  distinguished  from  the  scattering  coefficient - 
That  hope  has  not  been  realized  primarily  because  scattering 
overwhelms  radiant  energy  absorption.  Aside  from  either  pre¬ 
paring  a  sample  free  from  scattering  centers,  or  by  analysis 
of  emissivity  (or  reflectivity)  measurements,  we  have  not 
achieved  a  simple  reliable  method  of  determinina  the  effective 
absorption  coefficient.  Since,  as  reported  in  Part  II1  '  the 
absorption  coefficient  is  a  function  of  both  wavelength  and 
temperature,  determination  of  isothermal  emissivity  (or 
reflectivity)  remains  the  most  effective  way  of  fixing  that 
variable.  Once  this  is  done  however,  we  now  have  a  sound 
basis  for  evaluating  the  effects  of  both  microstructure  and 
temperature  gradients  on  the  radiant  energy  emission  over  a 
range  of  conditions  tested  by  both  analytical  calculations 
and  experimental  measurements. 

As  a  final  underpinning  to  confident  application  of  this 
approach  we  have  compared  the  Hamaker  approximation  and  an 
exact  solution  for  radiant  energy  emission  for  a  range  of 
thermal  conditions  and  real  material  constants  in  Section  6.0. 

We  have  found  that  the  Hamaker  approximation  is  suitable 
(within  0.1%)  for  the  conditions  tested. 
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APPENDIX  I 


1-  Characterization  of  the  Material  Specimens  used  for 

Transmissivity  and  Eteissivity  Measurements 

In  order  to  calculate  the  emissivity  and  transmissivity 
of  a  material,  certain  material  properties  must  be  character¬ 
ized-  Thes^  properties  fall  into  two  types:  (1)  properties 
which  are  unique  to  the  material  and  (2)  properties  which  are 
unique  to  the  sample*  The  index  of  refraction,  n,  is  an 
example  of  a  material  property,  while  the  scattering  coeffi¬ 
cient  S,  is  a  property  which  is  unique  for  each  sample. 

1.1  ai2o3 

Emissivity  and  transmissivity  measurements  were  made 
on  four  alumina  samples,  characterized  as  follows: 

1.1.1  AD-995 

Volume  Fraction  Porosity  (P)*  0-044 
Pore  Radius  (r)*  2.31*1 

1  Index  of  Refraction  Absorption  Scattering 

(n) *  Coefficient (cm- Coefficient (cm-1) 


a* 

a*=  2a 

S* 

1 

1-75 

0.38 

0.76 

385 

2 

1.74 

0-12 

0.24 

372 

3 

1.72 

o 

r-4 

• 

O 

0.20 

357 

4 

1-65 

0.24 

0.48 

327 

5 

1  .60 

2.6 

5.2 

277 

*These  symbols  are  the  same  as  used  throughout  the  text. 

Figures  1.1  and  1.2  are  photomicrographs  of  AD- 
995  alumina  before  and  after  testing. 

1.1.2  AD-85 

Volume  Fraction  Porosity  (P)  0.103 
Pore  Radius  (r)  2.31jj 


X 

n 

a 

S 

1 

1.75 

0.38 

348 

2 

1.74 

0.12 

348 

3 

1.72 

0.10 

331 

4 

1.65 

0.24 

278 

5 

1.60 

2.6 

178 

Figures  1.3  and  1.4  are  photomicrographs  of  AD- 
85  before  and  after  emissivity  measurements  were  run. 

1.1.3  AJ-1 

Volume  Fraction  Porosity  (P)  0.119 


Pore  Radius 

(r)  2 

X  n 

tt(cnT^) 

S  (cm--*-) 

1  1.75 

0.38 

973 

2  1.74 

0.12 

943 

3  1.72 

0.10 

907 

4  1.65 

0.24 

813 

5  1. 60 

2.6 

734 

Figures  1.5  and  1.6  are  photomicrographs  of  Ai-1 
taken  before  and  after  emissivity  measurements  were  run- 

1.1.4  Ai-4 


Volume  Fraction  Porosity  (P)  C.145 


Pore  Radius 

(r) 

X 

X  n 

e  (cm--*-) 

S  ( cm'^  ) 

1  1.75 

0.38 

2373 

2  1.74 

0-12 

2133 

3  1.72 

0.10 

1787 

4  1.65 

0.24 

1200 

5  1.60 

2.6 

787 

Figures  1.7  and  1.8  are  photomicrographs  of  the 
sample  before  and  after  emissivity  testing. 

1.2  SrTi03 

Two  specimens  of  strontium  titanate  with  different 
densities  were  fabricated  by  American  Lava  Company.  After 
thermally  annealing,  these  were  ground  to  the  required  shapes 
for  emissivity  and  transmissivity  measurements. 

1.2.1  SrTiO^-A  and  SrTiO^-B  A  B 

Volume  Fraction  Porosity  (P)  0.069  C.043 

Pore  Radius  (r)  1.6xp  2.8|i 
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A  B 


X 

n 

otcra'^l 

S  (cra“l) 

S  (cm“l) 

1 

2.31 

0.56 

1004 

549 

2 

7 

0.90 

1004 

549 

3 

2.23 

0.32 

1004 

549 

4 

2.18 

0.29 

803 

513 

5 

2.11 

2.4 

663 

495 

Figures  1.3  -  1.12  are  before  and  after  photomicro¬ 
graphs  of  the  strontium  t\ tana te  samples. 

1 . 3  Pyroceram 

Samples  of  Pyrcceram  9606  and  Pyroceram  9606  for 
emissivity  and  transmissivity  measurements  were  supplied  by 
Corning  Glass  Works.  The  9606  and  9606  Pvrocerams  were  both 
fully  dense  (i.e.,  no  porosity),  and  of  optically  indeter¬ 
minate  grain  size  so  that  no  characterization  of  them  could 
be  made  fcv  the  parameters  used  for  the  alumina  and  strontium 
titanate  specimens.  Figures  1.13  and  1.14  are  photomicro¬ 
graphs  of  Pyroceram  9606  and  9606,  but  there  is  no  detail 
visible. 
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Figure  1-1  Photomicrograph  of  AD-995  Before  Non-Isothermal 
Test  (400x) - 


Figure  1.2  Photomicrograph  of  AD-995  After  Non-Isothermal 
Test  (400x) . 
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Figure  1 . 3 


Photomicrograph  of  AD-85  Before  Non-Iso thermal 
Test  (400x) . 


Figure  1.4  Photomicrograph  of  AD-85 
Test  (400x) . 
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After  Non-Isothermal 


Figure  1.5  Photomicrograph  of  AX-1  Before  Non-Isothermal 
Test  (400x) . 


Figure  1.6  Photomicrograph  of  AX-1  After  Non-Isothermal 
Test  (400x) . 
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1.8  Photpmicrograph  of  AjJ-4  After  Non -I  so  thermal 
Test  (400x)  . 


-87- 


Figure  1 . 9 


Photomicrograph  of  SrTiC^A  Before  Non-Isothermal 
Test  (400x) . 


Figure  I . 10  Photomicrograph  of  SrTiC^A  After  Non-Isothermal 
Test  (400x) . 
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Figure  I. 11  Photomicrograph  of  SrTiC^B  Before  Non-Iso thermal 
Test  (400x) . 
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Figure  1.12  Photomicrograph  of  SrTiC^B  After  Non-Iso thermal 
Test  (400x) . 
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APPENDIX  II 


Spectral  Emissivity  of  Radiating  Energy  from  a  Radiating 
Material  with  an  Unknown  Surface  Temperature  and  a  Known 
Gradient  in  the  Body 

We  shall  assume  that  a  sample  of  radiating  material  has  a 
temperature  distribution  given  by 


T(x)  =  T  +  (-—)  x  (  1) 

o  dx  o 

in  the  interior,  i.e.  sufficiently  far  from  the  boundary,  x 
measures  the  distance  into  the  material  from  the  boundary.  We 
wish  to  determine  the  temperature  distribution  near  x  =  0  and 
also  ascertain  that  Eq.  1  is  a  valid  representation  of  T  in  the 
interior.  For  now  assume  that  the  absorption  coefficient  a, 
the  scattering  coefficient  S,  and  the  index  of  refraction  n  are 
independent  of  frequency.  We  have  the  following  equations  for 
the  flux  ttF(x)  and  the  mean  intensity  I*  (x)  . 

■dJixl  =  _  4a  f  (X)  +  4<*n"a  T4  (x)  ;  (  2) 


and 


4  di  (x) 
3  dx 


(a  +  S)  F  (x)  . 


(  3) 


These  equations  may  be  combined  to  obtain 


iLl-T-l*!  -  3a  (o'  +  S)  I  (x)  =  -3a  (a  +  s)n2a  T4(x)  .  (  4) 

dx 


Now  the  sum  of  the  energy  flux  due  to  conduction  and  the 
energy  flux  due  to  radiation  must  be  a  constant.  Therefore 


k 

dx 


+  rrF(x) 


=  H, 


(  5) 
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where  k  is  the  thermal  conductivity  of  the  material  and  H  is 
the  total  energy  flux.  Mote  that  H  <  0  for  energy  coming  out 
of  the  material.  Substituting  F(x)  from  equation  3  into 
equation  5  we  find 


_k  dT(x) 


3(or  +  s)  dx 


H 


(6  ) 


where  A  is  a  constant.  Therefore, 


I  (x)  =  *  3"^4^  ~S^  [hx  +  kT(x)  +  Aj.  (8  ) 


Substituting  from  equation  8  in  equation  4,  we  obtain 


rd  T  1  2  ‘  4 

k| — ~  -  3a  {a  +  s)  Tj  =  -4TTan  a  T 

'"dx 

3a (a  +  s) (Hx  +  A) .  (9  ) 


First  check  that  equation  1  is  a  valid  interior  solution  for 
the  temperature  field.  Let 

4  4  3  dT 

T  =  T  +  4T  (—4  -x  (10) 

o  o  dx  o 

ia  the  interior,  consistent  with  equation  1.  Then  equation  9 
becomes 

-3a  (a  +  s)  kf  T  +  (”)  x]  = 

L  o  'dx  o  J 

2  '  r  4  3  dT  1 

-4TT«n  a  I  T*  +  4T  (~)  x 
L  o  o  'dx  c  J 

+  3a (a  +  S)  (hx  +  A)  .  (ll) 
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Therefore 


H  =  -k(^-) 
dx  o 


16r 


3(»  +  S) 


2 

n  a 


_3  ,dT. 

T  (— ) 
o  dx  o 


(12) 


and 


A 


4tT 

VTTl  J _ - 

o  3(»  +  S) 


2  '  m4 

n  c  T 

o 


(13) 


Hence  the  observation  of  T  at  two  interior  points  of  the 
material  determine  A  and  H,  since  two  observations  give  Tq  and 

/dT 

diT~)  o *  We  ^ave  verifie<^  that  equation  (1)  is  a  valid  interior 

solution  for  T.  Now  we  wish  to  find  out  how  T  varies  near 
x  =  0.  Let  the  total  temperature  field  T  be  written  as  the  sum 
of  an  interior  and  a  boundary  contribution.  Write 


T 


* 

T 


+  T 

o 


+ 


(f)  X 

dx  o 


(14) 


Substituting  equation  (14)  for  T  in  equation  (9) ,  and  using 
equations  (12)  and  (13)  for  H  and  A,  we  obtain  the  following 
equation  for  T* . 


2  '  3 

,2  *  16TTcm  a  T 

d  T  r  o  -»  * 

- - —  -  I  3®(Q'  +  S)  -  - - - T  =  0  (15) 

,  -c  L.  K  J 

dx 


We  have  neglected  terms  which  are  quadratic  in  the 
*  ,dT 

quantities  T  and  {- — )  ,  which  are  assumed  to  be  small.  The 

dx  o 

validity  of  this  approximation  must  be  verified  a  posteriori . 

2  16nQ'n2a  ’  T2 

If  L  =  3o'{o'  +  s)  -  - 7 -  is  positive,  then  the  solution 

*  -Lx 

of  equation  (15)  which  approaches  0  as  x — *  ®  is  T  =  Ce  ,  where 
C  is  any  constant.  So  the  temperature  field  is  given  by 

T(x)  =  T  +  (~)  x  +  Ce  ^  (16) 

o  dx  o 
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Using  equation  8,  we  find 


i  <X)  “  3t"*n  S)  K  +  A  +  K  [To  +  (S>oX  +  Ce’LX]}, 


(17) 


and  from  equations  (3)  and  (17)  we  find 


p'txi  "HH  +  K[<S>c  -  <^_LX]}  (18) 

Now  C  may  be  evaluated  by  imposing  the  boundary  condition 

l'+(o)  =  p  I~(o) #  (19) 

where  p  is  a  mean  reflection  coefficient.  Now 

l'  +  (x)  =  I  (x)  +  (20) 


and 

I  (x)  =  I  (x)  -  . 


So  (19)  is 


I  (o)  + 


F(o) 

2 


P  (I  (o) 


IlsL 

2  ’ 


or  (1  -  p)l'(o)  +  (“-^(0)  =  0 
Using  (17)  and  (18) ,  (22)  becomes 


(21) 


(22) 


a  -  p>  lia-4vJl  {a  +  K  [T0  +  C]} 

H*  H"  +  4'i>o  -  -]}  -  0 


+ 


1 
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This  equation  gives 


C 


(7 - f)  (71)  + 

-L  y,  v^X  O 


(f 


1  t_£ 


3  [O'  +  S)  i  -  C 


H 

K  J 
L 


1 


(23) 


Now  we  use  this  temperature  distribution  (16)  tc  in¬ 
vestigate  the  radiation  field  as  a  function  of  frequency. 
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We  shall  assume  that  a  sample  of  radiating  material  has 
a  specified  temperature  distribution  given  by 

d'T  _ T  - - 

T(x)  =  T  +  (-—)  x  +  Ce  (16) 

o  dx 

o 

where  C  is  given  by  equation  (23) . 

Let  the  sample  have  absorption  coefficient  av ,  scattering 

coefficient  S^,  and  index  of  refraction  nv .  These  may  all  be 

frequency  dependent.  Let  I^(x,  9)  be  the  specific  intensity 

of  radiation  at  frequency  v  in  the  material,  where  8  is  the 

polar  angle  with  respect  to  the  +x  axis.  Assume  axial 

symmetry.  The  equation  of  transfer  for  the  radiation  in  the 

material  is  dl'(x,  8)  , 

cos  9  — ~-r* -  =  -(«  +  S  )I  (x,  9)  +  j  (x)  .  (24) 

dx  v  v  v  v 

j^x)  is  the  emission  coefficient  for  radiation  at  frequency 

v  in  the  material.  j^(x)  has  two  terms,  one  due  to  thermal 

emission  at  frequency  v  and  the  other  due  to  scattering  of 

radiation  at  frequency  v.  j^,(x)  is  given  by 

2  S  rr  , 

j^fx)  a  Q'i^nyBi/(T(x) )  +  -|  f  I^tx,  9)  sin9d9  (25) 

w  o 


Bj,(T)  is  the  Planck  function  given  by 


B  (T)  = 
V  v 


1 

hi'/KT 


(26) 


~  -1 


where  h  =  Planck's  constant,  £  =  Boltzmann's  constant,  and  c 
is  the  velocity  of  light  in  vacuo. 

dT  —Lx 

In  accordance  with  T(x)  =  T  +  (~)  x  +  Ce  we  write 

o  dx  c 

VT<X>>  =  Clo>  +  (^[(£-)0X  +  Ce_IjX]  (27) 

o 


It  is  convenient  to  replace  the  variable  9  by  p  =  cos9 
in  the  following  analysis.  So  we  will  now  write  I^(x,  p)  to 
denote  the  specific  intensity  of  radiation  at  frequency  v,  at 
a  depth  x  in  the  .material,  in  the  direction  9  =  cos“l  p.  The 
equation  of  transfer  (24) may  now  be  written  in  the  form 
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(28) 


dl^x,  p) 
dx 


-(or  +  S  )I  (x,  M.)  +  a  n  B  (T(x))  + 
'  V  V  V  v  V  V 


I^(x,  p)dp 


From  now  on  we  will  suppress  the  subscript  v  on  &v,  Sv, 
nv,  and  B^fT) .  But  it  is  important  to  remember  that  whenever 
a,  s,  n,  and  B(T)  appear  they  refer  to  these  quantities  at  the 
particular  frequency  v  with  which  we  are  dealing.  We  now  solve 
the  equation  of  transfer  (28)  by  the  method  of  discrete 
ordinates  [see  Chandrasekhar,  R.  T. ,  Chapters  II  and  III]  ^  . 
Divide  t^e  radiation  field  into  2m  streams  in  the  directions 
P-  (i  =  —  1,  —  2,  ...^m  and  ^  =  -p_^)  .  Let  any  integral  of 

the  form  ■  "*’  f(u)dp  be  aoproximated  by 
J-1 


m 

|  f  (M)dP  -  E  a  f(P  )  . 

J  J  J 

j  =  -  m 

j  /  o 

We  will  use  a  Newton-Cotes  quadrature  formula,  with 
j  =  l,...m.  The  a.  are  chosen  so  that  the  formula 


(29) 

2.-1 


m 

f(H)dM  -  I  a_:  f  (P  ■ ) 
-1  j  =  -m  J 

j  /  o 


(30) 


is  exact  whenever  f(P)  is  a  polynomial 

In  qeneral  we  will  have  a.  =  -a  .  ,  and 

3~3 


of  degree  less  than  m. 


m 


1. 


For  m  =  3,  M  =  1/5,  P2  =  3/5,  ^  =  1,  and  a1  =  50/144, 

75/144,  and  a 3  =  19/144.  (See  Abramowitz  and  Stegun  ^  , 
886,  formula  25.4.15) .  For  m  =  4,  =  1/7,  ^  =  3/7 , 

5/7,  p  =  l,  and  a  =  0.34595,  a  =  0.15313,  a  =  0.41400, 

4  1  Z  (9) 

0.08692  (See  Kourganoff,  page  112)'  .  Writing  the 


equation  of  transfer  for  the  direction  p^  we  have 

dl  (x,  )  ,  2  g  ®  i 

-  *  -  (or  +  Sjl^x,  P.)  +  Q'n  B (T (x) )  +  -  Ea^l^x,  Mj) 

1  j  =  -m 

j  t  o 

i  -±1.  ±  (31) 

These  are  2m  first  order  linear  differential  equations  for  the 
2m  quantities  I^(x,  ^i) ,  i  «  +  1,  ...  +  m. 


Define  B°  =  B(T  )  and  B  -  (|f)  T 

O  pT  O 


Using  the  shortened  notation,  I^(x,  P^)  =  1^,  these 
equations  may  be  written  as 


VS  =  -  S)I.  +  an^{ B°  +  B'[,f,ox  +  Ce-Lx]] 


T  ,dT. 


(33) 


s  m 

2  E  ajIj* 

j  =  -  m 
j  /  o 

A  particular  solution  is 

i.  =  n2{[B°  +  B'  (f)0(x  -  ?rr)]  +  ffB'cv_Lx}  a  < 

S  in  j  j 

where  h^  =  D/a  +  S  -  LP  j  and  D  =  l/[  1  -  ~  I  '+  s  “  Lp  -  ** 

J  ^  j=-m  3 

To  this  must  be  added  the  general  solution  of  the  homogeneous 
system 

dl .  m 

p,-~  =  -  (or  +  S)  I .  +  f  I  a I  .  ( 


Vta  “  -  +  S;ii  +  2  L 

3=-m 

j^o 

We  seek  solutions  of  the  form 
+K  (O'  +  S)  x 

I.=g.e  p  i  =  +  1 ,  . . . +  m . 

l  i  —  — 


(35) 


(36) 


Substituting  this  expression  (36)  for  1^  into  the  differ¬ 
ential  equation  (35) ,  we  find 
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g.  (*  +  s)  (1  +  K  u.) 

i  pi 


s  m 

2  ~  3j9j 

3  =  -  m 

j  r  o 


(37) 


The  right  hand  side  of  this  equation  is  a  constant  in¬ 
dependent  of  i.  Therefore 


9i  1  +  K  P. 

P  i 


•  i  =  +1 1  '  •  *  * +m , 


(38) 


where  A  is  some  constant.  Substituting  this  back  into  the 
equation  (37)  for  g^  and  dividing  by  A  we  obtain 


a  +  s  = 


2  1  +  K  p. 

P  3 

3  =  -  m 
j  r  o 


(38) 


This  may  be  written  as 


1  =  -r 


1  S 


2  +  S 


m 

I 


a  . 
_JL 


1  +  d.K 
3  P 


3  =  -  m 

j  o 


recalling  that  =  "Pi'  a-i  ~  th^-s  becomes 


(40) 


1  = 


m 

I 


O'  +  S  .  ,  v2  2 

3=1  1  -  K  U 


(43) 


This  is  the  characteristic  equation  which  must  be  solved  for 
the  constant  K  .  If  S 


O'  +  S 


<  1,  the  equation 


1  = 


a  +  s 


m 

I 


a  . 

i 


2  2 
1  -  P.K 
3  P 

i  =  1 


(42) 
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has  2m  non-zsro  distinct  roots  which  occur  in  pairs  + 

(p  =1,  •  ••m) .  Suppressing  the  solutions  which  would 
let  1^  become  infinite  exponentially  as  x->  »  we  write 

-K  (o'  +  S) x 

p 

o  m  L  e 

1 - }  +  5 —  +  r  — E - 

+  s  b'  p£l  1  -  Kp  (43) 

+  acl.e'1*} 

1  J 

The  m  constants  L  are  determined  by  the  boundary  conditions 
at  x  =  0.  The  boundary  conditions  »t  x  =  0  are 


1^(0,  +  n)  «  p  (6)  1^(0, -ii)  ,  for  J 1  -  ~  <  v  <  1 

n 


(44) 


and 


1^(0,  +  u)  =  1^(0,  -  H)  for  0  <  (i  <  _  1_ 


n 


p  { 6 )  is  given  by 


(45) 


Pie) 


tan  (8  -  9  )  1 

2  '  J 

sin  (9  +  9  )  tan  (6  +  0  ) 


(46) 


where  |i=  cos  0  and  sin  8 


n  sin  8 


We  now  restrict  the  discussion  to  the  case  for  which  the 
four-point  (m  =  4)  Newton-Cotes  quadrature  formula  is  used  to 

r1 

evaluate  integrals  of  the  form  J  f(^)dp.  in  terms  of  the 

-1 

discrete  angles  the  boundary  conditions  (44,45)  are 


Iv(0,  Pi)  =  p(9i)I^(0,-^i)  for  ^  _  i. 


n 


<  <  1»  (47) 


and 
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(48) 


i’ (0,Mi)  =  1^(0, -^i)  for  0  <  < 


1 


2 

n 


For  the  Newton-Cotes  four-point  (m  =  4)  quadrature  formula, 
we  have  u  =  1/7,  s  3/7,  H  «  5/7,  M4  *  1.  If  n  >  1.43, 
we  have 

0  <  <  u2  <  M3  <  J  l  .  1_  - 

n 


Therefore 


I„(0,  n.)  =  I„(0,  -  Mi)  i  =  1,  2,  3,  (49) 


and 


V0'  u4>  =  P(0)I„(0.  -  ^4)  • 


(50) 


This  may  be  written 


Note  that 


I„(0,1)  =  p (0) 1^(0, -1) . 

.  2 


(51) 


p  ( 0)  ■  lim  o(9)-  lim-|  rsiMG-e  +  tan  (8  ,  -  9  J 
9^0  9-^0  "sin  {0  +  0  }  tan  (6  +  0  ) 


=  lim 


r,e  -  9  ,2i  , .  r 

( - r)  =  lim 


e  +  9 


e2U-V, 

n  1 

2  .  _  1.  2  J 

9-N)  6  (1  +  n* 


<*-=-£> 


(52) 


If  1.43  >  n  >  1.11,  0  <  p  <  ^  1  and 

n^ 


'  1  -  ^2  <  ^3  <  ^4  =  1 

n 
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In  this  case  the  boundary  conditions  are 


Iv  (0,  nt)  =  lv  (0,  -  ni)  i  =  1,  2 


(53) 


Iv  (0,  |J3)  «  p(03)Il/  (0,  -  \i3) 


and  1^(0,  M4)  *  p  (0)  1^  (0,  -  . 

Now  from  equation  (43)  , 


(55) 


o  m  L  e 


-K  {»  +  S)X 
P 


(x,  n.)  =  nV{(x  -  (~)  +  ~  +  £  - 

l  l  l  <*  +  S  dx  o  1  ,  1  - 

B  p=l 


P.K 
i  P 


.  O' CD  -Lx'i 

«  +  S-  LM.  6  J 
3 


:i(o.  V  =„VfT- 


,dT, 

^iWo 
a  +  s 


+  £  - - K - 

,  1  -  P.K 
p=l  l  p 


Q'CD  ^ 

O'  +  3  -  LP.  J 
"1 


(56) 


We  assume  for  now  that  n  >  1.43.  Then 


I.(C  p.)  =  1.(0,  P±)  i  =  1,  2,3 


i  (0,  i)  =  (r-T-~) 2  I.  (o,  -  l) . 
i  n  +  i  i 


Substituting  from  our  expression  (57)  for  1.(0,  p.)  we  have 

r  x 
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dT 

2  '  rB°  ^i'dx'o  “  Lp 

n  B  —  •  ~I  pSl  1  -  n.r  +  a 

3  .  .  .dT.  +  P 


^CD 


+  S  -  Lp 

1  p  1 


]  = 


O' CD 


2_tb°  .  Mi  fe’o  ;  , 

n  3  L  ‘  a  +  s  S1  +  U.K  +  a  +  s  +  LIJ.J 

B  p=l  l  p  1 

for  i  =  1,  2,  3 
and 


2  ’  rS’ 

n  B  !  ~T  - 


,dT. 

lT“)  in  L 


dx  o 


+  Z 


O' CD _ 1  _ 

a  +  S  -L  „ 


B 


a  +  S  “  1  -  K 

P®1  P 


(59) 


dT 

.  2  _  ,  o  (— )  m  L  _ 

(S-^4)  n2B  T“  +  -^E-2  +  s  - E_  +  _5SD -  1  (60) 

n  +  1  L1  <*  +  S  1  +  K  Q'  +  S+  LJ 

B  p=l  p 


These  equations  may  be  manipulated  to  obtain 


dT 

™  Lp  Kp _ ^dx^  o  _  a  CDL _ 

p=l  1  -  p2  K2  *  +  S  {<*  +  S)  2  -  lV 

i  p  i 

i  =  1,  2,  3 


(61) 


and 


m 

I  L 
p=l  1 


,  ,  .n  +  1.  2 

L_  .  -1—1  =  (ffl,  n  +  1 


Ll-K  'n  +  l;  1  +  KJ  dx  o  a  +  S 
P  P 

fn  -  1.  2 

.  ,n  -  1.  2~]  B  ^  ^  r  ‘n  +  l' 
vn  +  1  J  la 


+  S  +  L  a  + 


1 -  1 

S  +  L  J 

(62) 


Of  course  in  this  case  m=  4  and  the  K  are  the  4  positive 
roots  of  the  equation  ^ 


1 


S 

a  +  S 


4 

E 

3=1 


a  . 
_ 1 


,  2  2 
1  -  |i.  K 
J  P 


(42) 
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We  have  p^  »  1/7,  | =  3/7,  p^  =  5/7,  p^  =  1,  and 
a  »  0.34595,  a2  =  0.15313,  a3  =  0.41400,  and  a4=  0.08692. 

S 

The  values  depend  on  the  ratio  + "" .  Having  chosen 

a  and  S  we  may  solve  for  the  K  .  Then  the  coefficients  in 
the  4  linear  equations  (61,62)^for  the  Lp  are  known,  and  we 
may  solve  for  Lp.  Then  we  may  compute 

dT 

'  _ .  2  •  r  dx^  o  p  Lp  ,  Q'CD _ 1 

lu{0,  1)  "  n  B  |ff  +  s  .  S14-k  a  +  S  +  LJ 

B  p=l  p 

(63) 

The  intensity  of  radiation  emerging  from  the  sample  in  a  narrow 
cone  normal  to  the  surface  is 


'emergent 


(54) 


dT 

.  _  4nB _  r  ^dx^  o  B°  ™  Lp  +  Q'CD  ~i 

^emergent  =  (n  +  ^  2  L*  +  s  B'  p=1  1  +  Kp  o  +  S  +  lJ 

(65) 

The  reason  the  analysis  cannot  be  carried  out  explicitly 
in  general  is  that  the  characteristic  equation  (42)  for  K  is 
a  polynomial  of  moderately  high  order.  For  our  four-point 

2 

method,  equation  (42)  is  a  quartic  in  .  For  each  choice  of 

g 

-  this  equation  must  be  solved  numerically.  In  this 

a  +  S 

connection  we  can  make  the  following  observations,  Equation  (42) 
may  be  writcen  in  the  form 


F  (K  ) 
P 


S 

ot  +  S 


2  2 

j=l  1  -  p  K 
J  P 


It  is  clear  the  F(K-)  =  F (-K  )  and  F(O)  =  1  -  ■  —  -  = 

F  P  Oi  +  b 

a  1 

— - -  >  O.  Also,  F(K  )  is  undefined  when  K  =  —  ,  j  =1,  2,3,4 

a  +  s  p  P  Mj 
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Sr(X/>) 


In  fact,,  lim 


K 


F(K  )  = 
I_  P 
d 

i 


-  00  ,  and  lim 


Finally  we  note  that 


dF(K  ) 
_ _ 


dK 


F(K  \ 

K  P 

P  Mj 


=  +  00  . 


4  2o\p2K 

r,  IIP 

1 - TT  0  <  0  for  K  >  0 

j-1  (1  -  A2)2  p 

D  P 


From  these  facts  we  conclude  that  F  (K  )  has  the  form 

1? 


Thus  we  see  F(Kp)  has  four  positive  zeros,  with 


0<K  <  1<K0<~<K_<~<K.  <7 

1  2  b  3  3  4 


(57) 


That  these  are  all  the  positive  roots  may  be  seen  by 
multiplying  equation  (66)  by 


2  2 

P(K  )=  II  (1  -  q.  K  )  . 
P  j=l  3  P 


(68) 


P(K  )  is  a  quartic  equation  in  K  . 

P  P 

gives 


Multiplying  (66)  by  P(K^) 


p<y 


s 

a  +  s 


4  a  .  4 

Z  2 — ?  11  t1 

j=l  1  -  p  K  i=l 


2  2 
p.K 
l  p 


(69) 


This  may  be  simplified  to  give 

4  4 

oty  -  p(y  -  1  ai 11  (1  -  ►*!  Ko>  = 0  (70) 

F  j=l  i=l  p 


44 

2  2 

Now  Zajll  (1-p.K  )  is  the  sum  of  four  cubic  polynomials 

j=l  i=l  1  P 

2  2 
in  K  .  So  Q(K  )  is  a  quartic  polynomial  in  K  ,  and  the  equa- 
P  P  P 

tion 

Q(K  )  =  0  (71) 


has  four  real  or  complex  roots  for  K  .  But  the  roots  of 

P 

Q(K  )  =  O  are  the  same  as  the  roots  of  F(K  )  ■  O,  and  we  have 
P  p 

shown  that  F (K  )  has  four  positive  roots  K  .  So  the  eiqht 
P  (  P 

roots  K  of  Q(K  }  =  0  are  just  +K  ,  where  K  are  the  positive 
p  p  ~  P  p 

roots  of  F(K  )  =0.  To  find  the  actual  roots  numerically 
P  1 
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equation  (66 )should  be  reduced  to  the  form  of  equation  (70)  . 
Then  any  of  the  standard  numerical  techniques  for  quartic 
equations  may  be  used.  (See,  for  example,  Abramowitz  and 
Stegun,  page  20,  Example  6.)  The  inequalities  (67)  will  be 
useful  in  any  iterative  procedure  used  to  solve  (70)-  Once 
the  K  are  known,  solution  of  the  system  of  equations  (61 

and  62)  for  is  quite  straight-forward.  It  requires  only 

the  inversion  of  one  4x4  matrix. 

„  o  1 

In  solving  (62)  it  is  useful  to  note  that  B  and  B  are 

given  by 


B  = 


2hv' 


2  hi'/KT 
C  ~  o 

e 


(72) 


-  1 


and 


hiv 


B 


B 


KT 


1  -  e 


-hi'/KT 
~  o 


(73) 
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APPENDIX  III 


Exact  Solution  for  the  Transmissivity 
of  an  Absorbing  and  Scattering  Material 

For  a  beam  of  parallel  light  of  unit  intensity  that  is 
incident  normally  on  the  surface  of  a  specimen,  the  equations 
for  the  flux  are: 


~  =  sJ  +  s'j  -  (a  +  s)I  +  f'l 
dx  p  p 


j  T  i  I 

~=sl+sl  -  (a  +  s ) J  +  f  J 
dx  p  p 


dl 


=  -  (a  +  s  +  f  )  I  =  -y  I 
dx  P  P 


dJ 


dx 


=  -(a  +  s  +  f  ) J  =  -g'j 

P  P 


Solutions  of  these  equations  are: 

I  =  I  oxp(-g'x) 

P  P° 


J  =  J  exp(g  x) 
P  po 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


I  =  A{1  -  8  )exp(c  x)  +  B(1  +  8  )exD(-g  x)  -  13  I  exp(-g'x) 
o  o  o*o  po  r  3 


-V  J  exp(g  x) 
po 


(7) 


J  =  A ( 1  +  8  )exp(a  x)  +  B(1  -  S  )exp(-a  x)  -  V  I  exp(-g'x) 
o  o  o  o  po 


-U  J  ex^,(g  x) 
po  * 


(8) 


where  A  and 
ditions  and 


B  are  const  ints  determined  by  the  boundary  con- 
U  and  V  are  as  defined  in  section  4.1. 
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The  boundary  conditions  are: 


x  =  o,  I  =  (1  -  p  )  +  J  p  . 

p  no  p  m 

I  =  Jp  . 


x  =  d,  J  =  I  p  . 

p  P  ni 

J  =  Ip  . 

1 

Applying  the  parallel  radiation  boundary  conditions, 

x  =  d,  J  exp(g'd)  =  p  I  exp(-g'd)  (9) 

po  ni  po 


x  =  o,  I  =  (1  -  p  )  +  Jp  . 

pc  no  ni 


(10) 


(1-p  ) 

no 


P°  1  ~  P^iexp(-2g‘d) 


(ID 


(1-p  )  p  . exp(-2g'd) 

no  *ni  ^ 


P°  1  ~  pni  exP(~29  d) 


(12) 


Applying  the  diffuse  radiation  boundary  conditions 

x  =  o,  A ( 1  -  8  )  +  B ( 1  +  8  )  -  u  I  -VJ 

o  o  po  po 

p.  (a(1  +  8  )  +  B ( 1  -  S  )  -  U  J  -VI  \  (13) 

1  l  o  o  po  poj 

x  =  d,  A  ( 1  +  8  )exp(a  d)  +  B(1  -  8  )exp(-cr  d)  -  V  I  exp(-g  d) 
o  o  o  o  po 

U  J^oexp(g'd)  =  p^A  (1-8^)  exp(aQd)  +  B (1+8q) exp(-gQd) 

-  U  Ipoexp(-g'd)  -  V  JpQexp(y ’d)  j-  (14) 
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or 

a[(1  -  p.)  -  (1  +  p.)8o]  +  b[(1  -  p.)  +  (1  +  p.)80]  = 


I  (U  -  o  V)  +  J  (V  -  o  .U) 
po  1  po  1 


(15a) 


A  exp(o  d)[ (1  -  p  )  +  (1  +  p.)3  1  +  B  exp(-a  d) f (1  -  p  )  + 

O  1  1  OJ  O  L.  1 

(1  +  o.)3  1  =  exp(-g'd)  (V  -  p.U)  +  J  exp(g'd) (U  -  p .V) (15b) 
i  oj  po  l  po  1 


and,  solving  by  determinants,  the  denominator 


=  2f  (1  -  p  . )  +  (1  +  p  . )  8  "1  sinh  j  d  +  4(1  -  0  . )  8  cosh  a  d 

L  1  1  oj  o  10  o 

(16) 


A  =  I  {[Jpc(V  '  °iU)  +  "po10  '  piv)]  •  [(1  -  pi>  -  (1  +  Po)Bo] 
exp(-aod)  -  [jpo(U  -  p^V)  exp(g'd)  +  I  (V  -  P±U)  exp(-g'd)  ]  * 

[(1  -  p.)  +  (1  +  p.)0o]}  (17) 

B  =  J  {[jpo(TJ  -  piV)exp(g'd)  +  I  (V  -  p±U)  exp(-g'd)  j  * 

[d  -  p.)  -  (l  +  Pi)8o]  -  [Jpo(v  -  p.U)  +  Ipo(U  -  p.V)]  • 

[(1  -  Pi)  +  (1  +  Pi)S0]  exp(aQd)j-  (18) 
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Substituting  for  A  and  B  in  equation  (7) 

I(x=d)  -4  {[Jpo(V  -  piU)  +  Ipo(0  -  piV)](4V  + 

[Jpo(U  "  P±v)exP(g,d)  +  I  (V  -  piU)exp(-g'd)]  • 

P {  ( 1  -82)  -  o.(l  +  fi2)}  2  sinh  a  d  +  4  p.B  cosh  a  d~!l 
L  o  i  o  J  o  o  o  JJ 


-  U  I  exp(-g  d)  -  V 
po 


J  exp(g'd) 
pc 


(19) 


The  transmitted  energies  are: 
for  parallel  radiation 


=  (1  -  p  .)I  , 

•nr  p(x  =  d) 


(20) 


and  for  diffuse  I.  =  (1  - 


Po)][(x  =  d) 


(21) 


Since  we  have  considered  the  incident  radiation  to  be  of  unit 
intensity,  the  transmissivities  for  parallel  and  diffuse 
radiation  are 


T 

P 


(22) 


and 


(23) 
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APPENDIX  IV:  Description  of  the  Three  Black-Body  Radiation 
Laws 

The  three  black-body  radiation  laws  and  their  differences 
are  as  follows: 


1)  The  Stefan-Boltzmann  Law  is  for  the  total  emissive  power 
of  a  black  body: 

w  =  aT4  (lv 


2)  Wien's  Law  is  for  the  monochromatic  emissive  power,  at 
wavelength  X ,  of  a  black -body: 

X  5exp(-C2/XT)  (2) 

This  law  was  found  to  be  incorrect;  and  errors  at  long  wave 
lengths  are  considerable. 


3)  Planck  then  came  along,  and  corrected  the  expression 
derived  by  Wien  to  its  correct  form  for  the  monochromatic 
emissive  power  of  a  black -body: 


C1 


exp(C2/XT)  -  1 


(3) 
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